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Abstract 


Periodic  surfaces  are  the  microwave  analogs  of  optical  diffraction  gratings.  While 
diffraction  gratings  serve  as  spectrum  analyzers,  periodic  surfaces  have  found  many 
practical  applications  as  frequency  selective  surfaces  which  transmit  or  reflect  a 
preselected  band  of  frequencies.  In  examining  the  electromagnetic  scattering  from  a 
periodic  surface,  one  generally  assumes  the  surface  is  planar  and  is  truly  periodic  which 
implies  the  surface  is  infinite  in  extent.  Practical  considerations  make  deviations  from  the 
ideal  unavoidable.  For  example,  infinite  surfaces  are  impossible  and  true  periodicity  may 
not  be  achieved  due  to  manufacturing  tolerances  and  imperfections.  The  effects  of  such 
deviations  on  the  performance  of  a  given  periodic  surface  thus  are  matters  of  practical  as 
well  as  theoretical  interest. 

The  purpose  of  this  work  is  to  characterize  the  degradation  due  to  deviations  from 
strict  periodicity,  on  electromagnetic  scattering  for  plane  wave  incidence  upon  a  periodic 
surface.  To  render  the  problem  tractable  and  still  yield  useful  results,  the  two-dimensional 
case  of  a  strip  grating  consisting  of  infinitely  thin  conducting  metallic  strips  with  parallel 
edges,  separated  by  gaps  or  slits,  was  chosen  to  initiate  investigations  into  the  effects  of 
the  variations.  The  variations  are  assumed  to  be  small  random  perturbations  in  the  widths 
of  the  strips  and  slits  as  compared  to  an  ideal  periodic  strip  grating. 

Prior  works  in  random  linear  antenna  arrays  treat  the  element  location  as  a  random 
point  with  distributions  such  as  Poisson  or  uniform.  This  treatment  is  not  suitable  for  use 
as  a  basis  for  this  study  because  a  distribution  specifying  the  distance  between  successive 
points  is  more  appropriate  in  specifying  the  random  errors  in  the  width  of  each  strip  and 
slit.  Using  this  specification,  the  average  power  factor  of  a  linear  array  of  point  sources  is 
calculated  and  is  equal  to  the  power  spectral  density  of  the  random  process  specifying  the 
inter  element  spacing.  The  variance  and  standard  deviation  of  the  average  power  factor 


XI 


are  also  calculated  and  as  the  number  of  points  approaches  infinity,  the  standard  deviation 
approaches  the  mean. 

Random  variations  in  geometry  preclude  exact  analytical  solutions  using 
techniques  such  as  the  Wiener-Hopf  method.  Hence,  formulations  for  the  far-field 
statistical  average  power  pattern  of  a  strip  grating  with  errors  in  the  widths  of  the  strips 
and  the  slits  are  presented  using  two  approximate  methods.  The  first  formulation  utilizes 
the  Bom  approximation  in  which  the  unknown  aperture  fields  are  replaced  by  the  incident 
fields.  The  second  formulation  utilizes  an  approximation  for  the  unknown  fields  which 
satisfies  the  edge  condition,  i.e.  has  singularities  at  the  edges  of  the  slits/  strips. 
Approximations  for  the  scattered  fields  are  first  derived  using  perfect  electrical  conducting 
(PEC)  surface  equivalence  for  a  TEz  polarized  plane  wave  incident  upon  a  strip  grating 
consisting  of  an  infinite  PEC  screen  cut  by  a  number  of  infinitely  long  slits  (infinite  in  the 
z-direction).  Babinet's  principle  is  then  used  to  obtain  approximations  for  the  fields 
obtained  when  a  TMz  polarized  plane  wave  is  incident  upon  the  complementary  grating 
formed  by  interchanging  the  slits  and  strips  of  the  original  grating,  i.e.  a  number  of 
infinitely  long  strips  in  free  space.  Expressions  for  the  average  power  pattern  are 
developed  in  terms  of  the  following  variables;  number  of  slits  or  strips,  desired  width-  and 
spacing-to-wavelength  ratios,  and  the  characteristic  functions  of  the  probability  density 
functions  of  the  width  and  spacing  errors.  Examples  are  presented  showing  the  effects  of 
the  above  variables  for  both  uniform  distributed  errors  and  errors  with  a  density  function 
based  upon  a  cosine  function.  The  variance  and  standard  deviation  for  both 
approximations  are  examined,  and  ,  like  the  power  factor  of  the  array  of  point  sources,  as 
the  number  of  points  approaches  infinity,  the  standard  deviation  approaches  the  mean. 
The  results  are  compared  to  and  agree  with  the  average  of  a  number  of  patterns  computed 
for  realizations  of  gratings  containing  actual  randomly  generated  errors  in  the  width  of 
each  strip  and  slit. 
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TOLERANCE  THEORY  OF  PERIODIC  SURFACES 


I.  Introduction 


Purpose 

The  purpose  of  this  research  effort  was  to  develop  analytical  methods  to 
characterize  the  degradation  in  electromagnetic  scattering  when  a  plane  wave  is  incident 
upon  a  periodic  surface  having  small  but  random  variations  in  geometry.  The 
characteristics  of  interest  were  the  reflection  and  transmission  properties  as  functions  of 
angle  of  incidence,  frequency,  and  polarization.  To  render  the  problem  tractable  and  yet 
yield  useful  results,  the  analysis  considered  a  two-dimensional  periodic  surface  composed 
of  perfect  electrical  conducting  (PEC)  strips  separated  by  gaps  or  slits,  i.e.  a  strip  grating. 
In  the  ideal  case,  all  strips  have  equal  widths  and  all  slits  have  equal  widths.  This  research 
investigated  the  effects  of  small  changes  in  the  widths  of  the  strips  and  slits  from  the  ideal 
case. 

Background 

Periodic  surfaces  are  the  microwave  analogs  of  optical  diffraction  gratings.  While 
diffraction  gratings  serve  as  spectrum  analyzers,  periodic  surfaces  have  found  many 
practical  applications  as  frequency  selective  surfaces  which  transmit  or  reflect  a 
preselected  band  of  frequencies.  In  formulating  the  problem  of  reflection  and  transmission 
of  a  periodic  surface  within  the  framework  of  Maxwell's  equations,  one  generally  assumes 
the  surface  is  planar  and  is  truly  periodic  which  implies  the  surface  is  infimte  in  extent. 
Practical  considerations  make  deviations  from  the  ideal  unavoidable.  For  example,  infinite 
surfaces  are  impossible  and  true  periodicity  may  not  be  achieved  due  to  manufacturing 
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tolerances  and  imperfections.  The  effects  of  such  deviations  on  the  performance  of  a 
given  periodic  surface  thus  become  matters  of  practical  as  well  as  theoretical  interest. 

Similar  problems  arise  in  antenna  array  theory  and  have  been  investigated  by 
several  authors,  notably  Lo  (14)  ,  Ruze  (23),  Skolnik  (25),  and  Steinberg  (26).  In 
antenna  array  analysis,  the  location  of  each  element  is  treated  as  random  and  the  problem 
can  be  solved  using  well  established  time  series  formulations  if  one  replaces  distance  by 
time. 

A  number  of  authors  have  investigated  different  mode  matching,  variational,  and 
Wiener-Hopf  solutions  to  the  problem  of  the  exact  periodic  strip  grating.  Several  of  these 
include  Agronovich,  Marchenko,  and  Shestopalov  (1),  Baldwin  and  Heins  (5),  Chen  (8), 
Daniele,  Gilli  and  Viterbo  (9),  Ishimaru  (11:189-194),  Kieburtz  and  Ishimaru  (13), 
Luk'yanov  (16),  Luneberg  and  Westpfahl  (17),  Miles  (18),  Vanblaricum  and  Mittra  (28), 
Weinstein  (29:267-281),  and  Wu  (30).  No  work  has  been  done,  however,  on  gratings 
with  random  variations  in  geometry. 

Overview 

The  purpose  of  this  work  was  to  conduct  a  thorough  and  detailed  investigation  of 
the  effects  of  deviations  from  strict  periodicity  on  the  scattering  from  periodic  surfaces. 
To  make  the  problem  tractable  and  at  the  same  time  reflect  prevaihng  manufacturing 
standards,  the  deviations  were  assumed  to  be  small  random  perturbations  in  an  otherwise 
periodic  system.  The  two-dimensional  case  of  a  strip  grating  consisting  of  an  infinitely 
thin,  planar  PEC  screen  cut  by  an  infinite  number  of  periodically  spaced,  infinitely  long 
slits  with  parallel  edges  formed  the  basis  for  this  research.  The  effects  were  characterized 
in  terms  of  the  expected  value  of  the  far-field  power  pattern  of  the  scattered  fields. 

Before  starting  on  the  strip  grating  problem,  two  chapters  are  included  to  provide 
important  background  information  and  theory.  First,  the  average  power  pattern  of  a  linear 
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array  of  point  sources  with  random  spacing  between  points  is  determined  in  Chapter  11. 
Associated  with  this  chapter  is  an  appendix  which  presents  the  calculation  of  the  variance 
of  an  array  with  an  infinite  number  of  point  sources.  The  methods  developed  in  this 
chapter  and  the  associated  appendix  for  dealing  with  random  spacings  in  the  array  serve  as 
the  basis  for  analyzing  the  random  errors  introduced  into  the  strip  grating  in  later  chapters. 
Second,  basic  electromagnetics  theory  used  in  the  formulation  of  the  scattered  fields  is 
presented  in  Chapter  III.  Here,  PEC  surface  equivalence  is  used  to  derive  a  general 
expression,  in  terms  of  the  unknown  equivalent  magnetic  surface  current,  for  the  fields 
produced  when  a  TEz  polarized  field  is  incident  upon  a  strip  grating  with  the  slits  and 
strips  infinite  in  the  z-direction.  Babinet's  principle  is  then  used  to  derive  an  expression  for 
the  fields,  in  terms  of  the  unknown  equivalent  electric  surface  current,  produced  with  TMz 
incidence  upon  the  complementary  grating.  The  complementary  grating  is  obtained  by 
interchanging  the  slits  and  strips  in  the  original  grating.  As  a  result  of  Babinet's  principle, 
the  equivalent  magnetic  surface  current  for  TEz  incidence  upon  a  strip  grating  is  the  same 
as  the  equivalent  electric  surface  current  for  TMz  incidence  upon  the  complementary 
grating.  Consequently,  the  same  integral  equation  is  used  to  represent  the  scattered  fields 
for  both  cases. 

The  next  two  chapters  utilize  two  different  approximations  for  the  unknown 
equivalent  surface  currents.  Chapter  IV  presents  the  results  of  utilizing  the  Bom 
approximation  in  which  the  unknown  fields  are  replaced  by  the  incident  fields  while 
Chapter  V  presents  the  results  of  utilizing  an  approximation  which  satisfies  the  edge 
condition.  For  both  methods,  expressions  for  the  far-field  average  power  pattern  are 
derived  for  both  gratings  with  no  errors  and  gratings  with  independent,  identically 
distributed  errors  in  width  (i.e.  the  length  of  each  slit  for  the  TEz  mode  or  the  length  of 
each  strip  for  the  TMz  mode)  and  spacing  (i.e.  the  length  of  each  strip  for  the  TEz  mode 
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or  the  length  of  each  slit  for  the  TMz  mode).  The  probability  density  functions  (PDF)  of 
the  errors  in  width  may  or  may  not  be  equal  to  the  PDF  of  the  errors  in  spacing. 

The  far-field  average  power  pattern  is  computed  and  presented  in  graph  format 
using  uniformly  distributed  PDFs  with  zero  mean  for  the  following  conditions: 

1)  Width  errors  alone,  spacing  errors  alone,  and  both  width  and  spacing  errors 

2)  Width  equal  to  spacing,  width  less  than  spacing,  width  greater  than  spacing 

3)  Maximum  error  values  (expressed  as  a  percentage  of  desired  width  or  spacing) 

of:  0  (i.e.  no  error),  5, 10, 15  and  20% 

Graphs  are  presented  for  gratings  with  a  finite  number  of  slits/strips  to  allow  comparison 
between  the  patterns  of  gratings  with  errors  to  gratings  with  no  errors.  In  addition,  graphs 
are  presented  for  gratings  with  an  infinite  number  of  slits/  strips.  Finally,  to  provide  a 
more  realistic  picmre  of  the  average  power  pattern,  patterns  are  provided  for  PDFs  based 
on  a  cosine  function  which  distributes  most  of  the  errors  around  the  mean  value  of  zero. 
Associated  with  these  two  chapters  are  appendices  which  present  the  calculation  of  the 
variance  of  a  strip  grating  with  an  infinite  number  of  slits/  strips  for  both  the  Bom 
approximation  and  the  edge  condition  approximation. 

To  validate  the  expressions  for  the  average  power  pattern  derived  in  Chapters  IV 
and  V,  the  results  of  a  number  of  trial  realizations  are  presented  in  Chapter  VI.  Here, 
errors  in  a  number  of  gratings  are  randomly  generated  using  the  uniform  or  cosine  density 
functions.  The  average  power  pattern  is  computed  for  each  grating  and  the  results 
averaged  together.  This  average  of  the  realizations  is  then  compared  to  the  statistical 
average  computed  using  the  expressions  derived  in  Chapters  IV  and  V .  Plots  showing  the 
standard  deviation  of  the  realizations  are  also  provided. 

Finally,  the  results  of  this  research  are  summarized  in  Chapter  VII  which  outlines 
the  major  accomplishments  and  provides  observations  and  conclusions  arising  as  a  result 
of  this  effort 
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II.  Linear  Array  of  Point  Sources 


Overview 


Before  looking  at  the  problem  of  the  strip  grating,  it  is  insightful  to  calculate  the 
average  power  pattern  of  a  hnear  aperiodic  array  of  point  sources.  This  is  equivalent  to 
the  average  power  pattern  of  a  linear  aperiodic  antenna  array.  As  stated  previously, 
several  authors  have  investigated  the  problem  of  random  errors  in  antenna  arrays.  In  the 
works  referenced  in  Chapter  I,  the  authors  calculate  the  average  power  pattern  of  a  linear 
aperiodic  array  where  the  location  of  each  element  is  treated  as  random.  The  distribution 
of  the  element  locations  is  typically  specified  by  normal  or  Poisson  PDFs.  These 
calculations,  however,  are  not  useful  in  this  investigation  because  they  consider  the 
location  of  each  element  as  the  random  quantity,  not  the  spacing  between  elements. 

Average  Power  Pattern  of  a  Linear  Array  of  Point  Sources 

Let  N+1  point  sources  be  arranged  on  a  line  with  the  location  of  the  n^^  source 
denoted  by  as  shown  in  Figure  2.1.  Define  the  distance  between  the  point  sources  to 
be  given  by: 


Ln=X„-X„. 


,n  =  l,2,...N 


The  array  factor  is  given  by  (15: 1 1-9): 

n=0 


(2-1) 


(2-2) 
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where 


u  =  sin  0 

0  =  angle  measured  with  respect  to  the  normal  to  the  array  (2-3) 

k  -  27t  /  X,  A,  •=  wavelength 

The  average  power  factor  will  be  defined  as: 


|f(u)|  = 


N+1 


N 

Ee 

n=0 


jkuX„ 


(2-4) 


For  the  case  where  the  spacing  is  periodic  with  period  T,  then: 


Ln  =  T  for  all  n 


(2-5) 


and  the  average  power  factor  becomes: 
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|2 


1 

N  +  1 


^  gjkunT 
n=0 


'>2 


(2-6) 


In  the  limit  as  the  number  of  point  sources  approaches  infinity,  the  average  power  factor 
becomes: 


lim  |fexact(u)|^  =  4  X 


N— >oo 


nX, 


(2-7) 


Now  let  the  array  have  random  errors  in  spacing  such  that  the  Lj,  are 
independent,  identically  distributed  random  variables  with  a  PDF  denoted  by  p(L).  The 
average  power  factor  is  now  given  by: 

i^  =  ^E{|f(u)f}  (2-8) 

where  E{|f(u)|^}  denotes  the  expected  value  of  |f(u)p given  by  (21: 138): 

E{if (u)f  }  =  |f (u)|V(L)dL  =  jf (u)pp(L)dL  (2-9) 

The  last  equality  arises  from  the  fact  that  p(L)  =  0  for  L  <  0,  i.e.  the  distance  between 
sources  is  always  a  positive  number.  Thus 
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|21 


1  N 

— y 

r  .  1  ^ 


JkuX„ 


1 


N  +  lnt^o 

■  N  N 


N  +  1 
1 

N  +  1 


jkuXn  -jkuX„. 


ln=0  m=0 


N  N-n 


EUN  +  l)  +  2Re£  X® 


jku(Xu+,^  X|j^) 


n=l  m=0 


(2-10) 


Where  "Re"  denotes  the  real  part.  But, 

X„™-X„=L„,+L„,+  -  +L„„  (nterms)  (2-11) 

and 


N  N-n 


I  n=l  m=0 


jku(Lm4.i+Lm+2‘*"  ”■  "^^m+n) 


O  N  N-n  ^  , 

=  !  +  £■{  “^-Re^  ^gjku(L„+i+L„+2+  -  +L„+„) 

^  n=l  ni=0 


N- 


(2-12) 


=  1  +  — ^ReV  V  ■"  ■^^"'+"^1 

N  +  1  SlTo 


Where  the  last  equality  arises  from  the  fact  that  the  expected  value  of  a  sum  of 
independent  random  variables  is  the  sum  of  the  expected  value  of  each  separate  randorn 
variable.  Since  all  the  Ljj  are  independent  but  identically  distributed, 


"^^m+n)  ^  —  E  |e  jkuL„+i  jg|gjkuL„+2  j. . .  E{el'^“^  m+n  ^ 

=  E"{el‘^“^} 


(2-13) 


and 
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(2-14) 


N  N-n 


^  ^  n=i  m=0 

=  1  +  -^ReV  (N  + 1  -  n)E"  {ej‘'“4 
N-f-1  ^ 


Now  define  the  characteristic  function  0(u)  as: 


0(u)  =  pp(L)ej'^“^dL  =  pp(L)eJ^dL  (2-15) 

J— oo  Jo 


The  characteristic  function  has  the  following  general  properties  (21:153-154): 


^+00 

1.  0(0)=  [  p(L)dL  =  l 

J^oo 

2.  |0(u)|  =  |j''V“^p(L)dL  <  J%(L)dL  =1  |0(u)|<l 

3.  |0(u)l<l  for  u  0  unless  L  takes  on  values  forming  an 

arithmetic  progression(i.e.  periodic) 


(2-16) 


Substituting  for  the  characteristic  function,  the  average  power  factor  becomes: 


N 

l  +  2Re^O"(u)- 

n-1 


2 

N-l-1 


N 

Re^nO"(u) 

n-1 


(2-17) 


Now  use  the  relations: 


1-0 


and 


o+no^^^-(n -1-1)0'^^* 

(ITJ? 


(2-18) 


(2-19) 
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to  obtain; 


(2-20) 


Application  of  L'Hopital's  rule  yields: 


|f(0)f  =N+1 


(2-21) 


For  the  case  of  an  infinite  number  of  point  sources,  i.e.  N  — > ,  use  the  following 
relations  for  |<E>|  <  1: 


and 


to  obtain: 


limZ*”  = 

n=l 


1-0 


N 


N-^oo 


o 

(i-®r 


-lim 

N“>oo 


l-H2Re 


J  I 
ll-0(u)r 


=  1  +  2RC' 


f  »(u)  1 

'l-4>(u)J 


(2-22) 


(2-23) 


(2-24) 


For  the  case  where  u  =  0 ,  application  of  a  Taylor  series  expansion  around  u  = 

0  can  be 

used  as  follows: 

Let 

u  ^  p.+  jo) 

(2-25) 

then  retaining  only  the  first  order  terms  of  a  Taylor  series  expansion  for  small  p  and  o), 

gjkuL  1  + jk(p,+ j(o)L  =  l-Ci)kL+  jkpL 

(2-26) 

and 

<I>(p+  j(0)  p(L)[l  -  (okL  +  jkpL]dL 

=  l-(OkL  +  jkpL 

(2-27) 

where 

L  =  E{l}  =  JJ)(L)LdL 

Substituting  Equation  (2-27)  into  Equation  (2-24)  yields: 

i-T - ^2  l-ll-wUL  +  jkpLf  1- (l-col^f +(kMLf 

f(p+j(0)  -  '  _  2-  7 ”-x2  /  -72- 

|l-(l-a)kL  +  jkpL)|  (cokLj  +(knL) 

2a)-kL((0^  +  p^) 

kE(co^-i-p^) 

(2-28) 

(2-29) 

or  retaining  only  the  first  order  terms; 
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(2-30) 


2(0 

kL((0^  +  |x^) 


Using  the  following  relation  (27: 166): 


- T-->2n5(|i.)  as  (o->0 

0)“^  +n 


(2-31) 


Equation  (2-24)  becomes: 

=  ^  for  ^  =  0  (2-32) 

kL  U 

which  is  the  same  as  for  an  array  with  exact  periodic  spacing. 

The  variance  and  standard  deviation  (o)  of  the  average  power  pattern  of  an  array 
with  an  infinite  number  of  point  sources  is  calculated  in  Appendix  A  with  the  following 
results: 


a  = 


for  u  0 
for  u  =  0 


(2-33) 


Comparing  Equation  (2-6)  to  Equation  (2-20),  we  see  that  the  average  power 
pattern  for  the  array  with  random  spacing  errors  is  related  to  the  average  power  pattern 
with  no  spacing  errors  via  terms  involving  the  characteristic  function  of  the  spacing 
errors.  Similar  expressions  will  be  developed  for  the  strip  grating. 
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III.  Electromagnetic  Theory 


Overview 

This  chapter  presents  basic  electromagnetics  theory  required  to  develop  the 
expressions  for  the  fields  scattered  by  a  strip  grating.  Expressions  for  the  scattered  fields 
are  derived  in  terms  of  integral  equations  involving  unknown  equivalent  magnetic  or 
electric  scattering  currents. 


Integral  Equations  for  an  Aperture  in  a  Conducting  Screen 

The  following  forms  of  the  free  space  Maxwell's  equations  expressed  in  Gaussian 
units  with  e^“*  time  convention  will  be  used  in  all  derivations: 


VxH  =  jkE 
VxE  =  -jkH 


(3-1) 


Consider  a  plane  wave  incident  upon  a  planar  PEC  screen  containing  one  or  more 
apertures  as  shown  in  Figure  3.1.  Let  A  denote  the  aperture(s)  and  M  denote  the  PEC 
surface.  An  equivalent  problem  can  be  obtained  by  wrapping  a  Huygen  surface  around  the 
screen  and  using  PEC  surface  equivalence  with  the  appropriate  equivalent  electric  and 
magnetic  surface  currents  (4:329-334).  The  equivalent  surface  currents  are  given  by 
(4:330): 


Jseq=nxH 

Msec,=Exn 


(3-2) 
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Figure  3.1  Plane  Wave  Incident  upon  PEC  Screen  with  Aperture 


where  n  is  the  outward  normal  from  the  Huygen  surface  and  E  and  H  are  the  total  fields 

on  the  surface  of  the  screen.  Since  the  screen  has  been  replaced  with  a  PEC  plane,  the 
Jseq  electric  currents  do  not  radiate.  Also,  =  0  on  the  original  PEC  surface,  leaving 

a  non-zero  Mg^^  magnetic  current  only  in  the  location  of  the  former  aperture  as  shown  in 
Figure  3.2.  Mggq,andMg^2  are  given  by  (4:330): 

Ms,qi=-nixEA=-yxEA 

Mseq2  =  -02  X  Ea  =  y  X  Ea  =  “Mseql 


where 


Ea  =  total  E  in  the  aperture 
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Figure  3.2  Equivalent  Problem 


The  vector  potential  F  is  given  by  (4:279): 


F(R)  =  J^Ms(R’)G<,(R,R')da' 

where 


Go(R,R') 


^-jk|R-Rl 

47r|R-R'| 


is  the  free  space  Green's  Function  and 


|R-R'  =  ■y/(x  -  x')^  +  (y  -  y ')^  +  (z  -  z')^ 


The  scattered  fields  due  to  the  vector  potential  are  given  by  (4:259) 


(3-7) 


E^=-VxF 

H*  =  -— (VxEJ  =  -ivx(-VxF)  =  -jkF— rV(V-F) 
jk  k  k 

The  total  fields  are  given  by,  using  image  theory  (7:97): 


'E‘+E^+F(2Mseq,) 

-1 

o 

A 

(3-8) 

.E‘(2Ms„,) 

A 

O 

'H'+H'+ff(2Ms„,) 

,y>0 

(3-9) 

,y<0 

where  E"(2Mseqi)and  H®(2Ms^,)  are  the  scattered  fields  due  to  2Mseqi  and 
E"(2Mseq2)  and  H®(2Mseq2)  are  the  scattered  fields  due  to  2Mseq2-  The  factor  of  two  is 

due  to  the  image  of  the  magnetic  current  in  the  presence  of  the  PEC  plane.  E'  and  H'  are 
the  incident  fields  while  E*^  and  are  the  fields  reflected  from  the  PEC  plane.  From 
image  theory,  the  components  of  the  reflected  fields  are  given  by  (3:158): 

E;,  (x,  y,  z)  =  -E|,  (x,  -y,  z)  (x,  y,  z)  =  H;,  (x,  -y,  z) 

Ey(x,y,z)  =  Ey(x,-y,z)  Hy(x,y,z)  =  -Hy(x,-y,z)  (3-10) 

E[  (x,  y,  z)  =  -E^(x,  -y ,  z)  (x,  y ,  z)  =  (x,  -y,  z) 

Boundary  Conditions  for  an  Aperture  in  a  Conducting  Screen 

Due  to  the  S5mimetry  of  the  problem  and  the  vector  potential,  the  following 
relations  apply  to  the  scattered  fields  (12:436),  (7:97),  (6:559): 
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E®(x,y,z)  =  E'(x,-y,z) 

H*(x,y,z)  =  -H*(x,-y,z) 

Ey(x,y,z)=  -Ey(x,-y,z) 

Hy(x,y,z)  =  Hy(x,-y,z) 

(3-11) 

E^(x,y,z)=  E*(x,-y,z) 

H",(x,y,z)  =  -H*(x,-y,z) 

On  the  PEC  surface,  the  tangential  E-field  components  must  be  zero.  Thus,  on 

the  PEC 

surface  M, 

Ex 

=  E,=0 

(3-12) 

In  the  aperture  A,  continuity  of  the  fields  requires  that 

e:,=e'+e;+e:i 

HL=H‘+H^+H^xi 

e;,  =  e;+e;+e;, 

(3-13) 

E:,  =  Ei-HE;+E:, 

Using  Equation  (3-10)  and  the  fact  that  x= 

•0  in  the  aperture.  Equation  (3-13)  becomes 

F®  —  F® 

^x2  ~ 

e;,=2e;+e;. 

TJS  _  TjS 

^y2  “  ^yl 

(3-14) 

E:a=E:. 

Hl2=2K+Ki 

Using  Equation  (3-1 1)  at  x=0,  we  see  that  in  the  aperture, 

K2=K 

K2=K  (3-15) 
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Equations  for  TEz  Mode  Incident  on  Screen  with  Infinite  Slits 

Now  consider  a  TEz  polarized  wave  incident  upon  a  conducting  screen  in  which 
the  aperture  or  apertures  are  infinite  in  the  z-direction.  For  the  geometry  shown  in  Figure 
3.3,  the  TEz  incident  wave  is  given  by: 

H*  =  z  =  z  (3-16) 

Since  the  incident  magnetic  field  has  only  a  z-directed  component  which  is 
independent  of  z  and  the  screen  has  no  variations  in  z,  the  scattered  and  total  magnetic 
fields  also  have  only  a  z-directed  component  which  is  independent  of  z. 


The  scattered  magnetic  fields  can  be  found  using  Equations  (3-4)  through  (3-7)  assuming 
Mg  has  only  a  z-directed  component  which  is  independent  of  z  (4:698-699).  Thus 


F(x,y)  =  z  Ms,(x')^-n - 2-^ - ^dx'dz'  (3-17) 

Ja  47tv(x-x'r+y''-Kz-z'r 
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Use  the  following  integral; 


to  obtain 


dz  —  j7iH[,^\ka) 


F(x,y)  =  -z| Ms,(x') H?)(kV(x-x')'+y')dx' 


(3-18) 


(3-19) 


where  LI  is  the  set  of  lines  on  the  x-axis  where  the  apertures  cut  the  z  =  0  plane. 
Substituting  Equation  (3-19)  into  Equations  (3-7)  and  (3-9)  yields 


where 


H  =  H,z 


(3-20) 


and 


H*  +H^  J^^Ms(x')  H®(kV(x-x')'+y')  dx' 

I  Ms(x0  H®(kV(x-x')Hy2)  dx' 


,y>0 

(3-21) 

,y<0 


Mg  —  Mggqjj  Mseq2z  ^Ax 

=  the  X  component  of  the  E  field  in  the  aperture 


(3-22) 


The  E-field  components  are  found  by  substituting  Equation  (3-20)  into  Equation  (3-1)  to 
obtain: 
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(3-23) 


E  --1^ 

jk  ay 

E 

^  jk  ax 
E,=0 

The  unknown  quantity  Mg  can  be  found  by  solving  the  following  integral  equation 
obtained  by  substituting  Equation  (3-21)  into  the  boundary  condition  for  the  tangential 
component  of  the  H-field  in  the  aperture  (Equation  (3-15)) : 

J^jMs(x')  H|f\klx-x'|)dx'  =  ■^H^(x,0)  (3-24) 

which  is  valid  for  points  in  the  slit  regions. 

Equations  for  TMz  Mode  Incident  on  Screen  with  Infinite  Slits 

Now  consider  a  TMz  polarized  wave  incident  upon  a  conducting  screen  in  which 
the  aperture  or  apertures  are  infinite  in  the  z-direction.  For  the  geometry  shown  in  Figure 
3.4,  the  TMz  incident  wave  is  given  by: 

E‘  =  E^  z  =  ^  z  (3-25) 

Since  the  incident  electric  field  has  only  a  z-directed  component  which  is 
independent  of  z  and  the  screen  has  no  variations  in  z,  the  scattered  and  total  electric  fields 
also  have  only  a  z-directed  component  which  is  independent  of  z.  The  scattered  electric 
fields  can  be  found  using  Equations  (3-4)  through  (3-7)  assuming  M^  has  only  a  x- 

directed  component  which  is  independent  of  z.  This  leads  to  a  solution  for  the  scattered 
electric  field  in  the  form  of: 
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Figure  3.4  Plane  Wave  Incident  upon  PEC  Screen  with  Slit,  TMz  Mode 


E,=z^  (3-26) 

dy 

An  alternative  solution  can  be  found  directly  from  the  TEz  mode  solution  using  Babinet's 
principle.  Babinet's  principle  is  as  follows  (3:166)  and  (6:560): 

1)  Let  Ej  =  F  and  Hj  =  G  be  incident  on  the  screen  from  the  +y  side  and  Ej  and  Hj  be 
the  resulting  fields  on  the  -y  side. 

2)  Let  Ej  =  -G  and  Hj  =  F  be  incident  on  the  complementary  screen  from  the  +y  side 
and  Ej  and  Hj  be  the  resulting  fields  on  the  -y  side.  The  complementary  screen  is  formed 
by  interchanging  the  PEC  surfaces  and  the  apertures. 


Then:  _  _  _ 

Ei+Hj  =F 

Hi-E2  =  G 


(3-27) 
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Babinet's  principle  can  be  used  as  follows; 


Let  El ,  E j ,  Hi ,  Hj  be  the  TMz  fields  and  E2 ,  E^ ,  Hj ,  be  the  TEz  fields.  Then, 


and 


Ei=E,iZ 

H2=H,2Z 


(3-28) 


E|  =  H‘=F,z 


(3-29) 


From  Equation  (3-27), 

^zi=Ki--\  Ms(xOH?\kV(x-x')^y')dx'  (3-30) 

where  L2  is  the  set  of  lines  on  the  x-axis  where  the  PEC  surfaces  cut  the  z-0  plane. 
Although  Mg  is  numerically  equal  to  the  magnetic  scattering  current  found  in  the  TEz 
case,  it  now  can  be  considered  to  represent  the  electric  current  density  induced  on  the 
surface  of  the  strips.  This  leads  to  introducing  an  electric  scattering  current  Jg  by  the 
following  substitution; 


Mg  — >  Jg 


(3-31) 


With  this  notation.  Equation  (3-30)  becomes 

Ezi  =  Ki  "T  f  Js(x')  H®(kV(x-x')^y')dx'  (3-32) 

2 


Thus,  for  the  TMz  mode. 
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E  =  E,z 


(3-33) 


where 


E,  =  E‘  H<„«(kV(x -X')" ♦y' )dx'  fory<0 

=  E‘+Ef°‘ 

The  E®’“‘  notation  is  used  to  differentiate  the  field  scattered  by  the  complementary  grating 
for  the  TMz  case  from  the  field  scattered  by  the  original  grating  for  the  TEz  case  as  given 
in  Equation  (3-8).  To  solve  for  the  electric  field  for  y>0,  superimpose  the  image  problem 
of  -E‘(x,-y,z)  incident  on  the  negative  side  of  the  screen  which  yields  -E^(x,-y,z) 

(3:163-164).  Now 

E2(x,y,z)-E^(x,-y,z)  =  E3(x,y,z)  +  E®’‘°‘(x,y,z)-E^(x,-y,z)-Ef°‘(x,-y,z)  ^2.35) 

=  E^(x,y,z)-E^(x,-y,z) 


since  E®’*°‘  is  an  even  function  in  y.  Thus  for  y>0, 

E^  (x,y,  z)  =  E^  (x,  -y,  z)  +  E^  (x,  y,  z)-  E^  (x,-y,  z) 

=  E^  (x,  -y,  z)  -  1^2  Js  (k-^(x-x')^  +y^  )dx' + E^  (x,  y,  z)  -  E^  (x,  -y ,  z) 

(3-36) 

Equations  (3-34)  and  (3-36)  have  the  same  form ,  thus 

E,  =  E‘  Js(x')  HS,^^(k-J(x-xO^+y^)dx'  for  all  y 

=  E‘-HEr 
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Note  the  similarity  in  the  form  of  the  integrals  for  the  TMz  case,  Equation  (3-37)  and  the 
TEz  case  Equation  (3-21).  The  H-field  components  are  found  by  substituting  Equation 
(3-33)  into  Equation  (3-1)  to  obtain: 


1  aE, 


jk  ay 


H 

y  jk  ax 
H,=0 


(3-38) 


The  unknown  quantity  Jg  can  be  found  by  solving  the  following  integral  equation  obtained 
by  substituting  Equation  (3-37)  into  the  boundary  condition  for  the  tangential  component 
of  the  E-field  on  the  (PEC)  surface  (Equation  (3-12) ) : 

Js  (xO  -  x'l)  dx'  =  ^  E^  (x,  0)  (3-39) 


which  is  vahd  for  points  on  the  original  PEC  surfaces.  Note  that  Equation  (3-39)  has  the 
same  form  has  Equation  (3-24)  obtained  for  the  TEz  case.  In  fact,  since 

E^  for  the  original  grating  =  for  complementary  grating 

L2  for  the  original  grating  =  for  the  complementary  grating 


then  as  expected  from  Babinet's  principle, 

Jg  for  the  original  grating  =  Mg  for  complementary  grating  (3-41) 


and  in  the  region  behind  the  grating,  i.e.  y<0. 
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for  the  original  grating  =  H*  for  complementary  grating 


(3-42) 


Equations  for  TEz  and  TMz  Modes  Incident  on  Strip  Grating 

The  integrals  for  the  TEz  mode  in  Equation  (3-21)  and  the  TMz  mode  in  Equation 
(3-37)  both  are  of  the  form: 


S  =  -  j  I,  H®(k7(x-x')^+y^)dx'  (3-43) 

where  S  is  the  scattered  field,  Ig  =  Mg  or  Jg,  and  L  =  Lj  or  L2  as  appropriate.  Thus  both 
can  be  evaluated  in  the  same  manner.  Now  consider  the  screen  to  contain  a  number  of 
strips  and  slits  infinite  in  the  z-direction  as  shown  in  Figure  3.5. 


For  the  TMz  case,  the  Wjj  represent  the  width  of  the  strips  and  the  a^  represent  the 
width  of  the  slits.  For  the  TEz  case,  the  an  represent  the  width  of  the  strips  and  the  Wq 
represent  the  width  of  the  slits.  For  future  discussion,  the  term  width  will  refer  to  the  Wn 
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while  the  term  spacing  will  refer  to  the  ajj.  Using  these  specifications  for  the  grating, 
Equation  (3-43)  becomes: 


S  =  r"l3(x') 


(3-44) 


where  n  ranges  from  Oto  +o°.  The  integral  equations  incorporating  the  boundary 
conditions,  Equation  (3-24)  and  Equation  (3-39),  become: 

y  i  r’l.(x')  H®(k|x  -  x'|)dx'  =  e'“"»  (3-45) 

^  2  •'x„ 

n 


Equation  (3-45)  is  valid  for  <  x  < 

Referring  to  Figure  3.5,  the  following  relations  hold  for  a  grating  with  no  errors: 

a„=A 

w„  =  W  (3-46) 

A-hW  =  T 

where  T  is  the  period  of  the  grating.  Setting  the  origin  at  Xq,  i.e.  Xq  =  0,  Equation  (3-44) 
becomes: 


S  =  H?\k7U^^^^fr7)dx' 

n  ^ 


(3-47) 


Using  the  change  of  variables 
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x"  =  x'-nT 


(3-48) 


Equation  (3-47)  becomes: 

S  =  X- Is  (x"  +  nT)  (k^(x-x"-nT)Sy^  )dx"  (3-49) 

n  2 

Since  the  grating  is  periodic,  the  current  in  the  slits  or  strips  is  related  according  to 
Floquet's  theorem  by  (24: 14); 


I,(x"-^nT)  =  l3(x")ej‘“'^*‘"® 


(3-50) 


This  leads  to 


S  =  Is(x") ll^\k^(x-x"-in'f+y^)dx"  (3-51) 

n  2 


The  equation  for  the  boundary  conditions,  Equation  (3-45)  becomes: 


Jj"  i^(x")  H^^)(k|x  -  X"  -  nT|)dx"  =  e^^ 


(3-52) 


For  the  purposes  of  this  analysis,  errors  in  geometry  of  an  ideal  grating  will  be  introduced 
as  perturbations  from  the  desired  values  of  each  strip  and  slit  width.  This  can  be 
expressed  as: 


subject  to  the  conditions: 


a„=A  +  e„ 

Wn=W  +  ^„ 


(3-53) 
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e{8}  =  e{|1=o 


(3-54) 


where  the  are  independent,  identically  distributed  random  variables  with  a  PDF  denoted 
by  pe(e).  Similarly,  the  are  independent,  identically  distributed  random  variables 

with  a  PDF  denoted  by  p^(^).  The  PDFs  pe(e)  and  p^(^)  may  or  may  not  be  equal  to  each 
other. 
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IV.  Born  Approach 


Overview 

This  chapter  presents  results  obtained  when  the  unknown  fields  in  the  integral  field 
equations  are  replaced  by  the  incident  fields.  This  is  an  application  of  the  Bom 
approximation  (20:1073)  and  is  equivalent  to  a  physical  optics  approximation.  It  is  valid 
in  the  high  fi’equency  region  where  the  wavelength  of  the  incident  wave  is  much  smaller 
than  the  widths  of  all  the  slits  and  strips.  Consequently,  it  is  especially  useful  in  problems 
of  optical  diffraction. 

Far  Field  Pattern  for  Grating  with  No  Errors 

For  a  grating  with  no  errors,  the  scattered  fields  are  given  by  (Equation  (3-47)): 

S  =  Sf ir'"  H?^(k70^^^^^f^)dx'  (4-1) 

n 

Using  the  relation  _ 

Hq^^x)  » asx-^oo  (4-2) 


the  scattered  fields  in  the  far  zone  become: 


Sic  fnT+W,  .  , _ 

-f  i,(x')  nr 

2Jm  ■J/«k 


^(x-x')^+y^ 


-/k^(x-x'f +y^-ii/4l 

*  ^  ^dx' 


(4-3) 


Next,  use  the  Fraunhofer  approximation: 
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R  for  amplitude  purposes 

R  +  x'  sin  0  for  phase  purposes 


(4-4) 


^(x-xO^+y^  «| 

where 

■Jx^+y^  =R 


(4-5) 


and  9  is  the  observation  angle  as  defined  in  Figure  4.1  to  obtain: 


(4-6) 


where  the  superscript  "f '  refers  to  the  far-field  approximation  for  the  scattered  fields. 


Figure  4.1  Incidence  and  Observation  Angles 


Consider  TEz  incidence  with  the  incident  electric  field  defined  by: 


jji  _  gjk(xsinei -t-ycosOi) 


(4-7) 


then. 
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I,(x)  =  M,  =-yxEy^^ 


-I  if^H, 

A-  I  >v  •*■  7  ^  7  ^ 

-vxE  =-yx - - y 

'y=®  jk  1  ^  5x 


z  =  cos0ie^‘‘’^“”®-z 

jk  ^  y=0 


I,(xO  =  coseie^‘“'“"®‘ 


for  both  the  TMz  and  TEz  cases.  Inserting  this  into  Equation  (4-6)  yields: 


s'  =  cos9i,/i  e f"*"'  e'^'dx' 

V2n  VR  „  " 

=  cos9,Ji//^  ^irr 

=  cos9.JX'e'%  e'™^ 

'vrit  Vr  u  y 


(4-10) 


where 


u  =  sin  0;  -  sin  0 


(4-11) 


Let  n  range  from  0  to  N,  then  Equation  (4-10)  becomes: 


Sf-cos9  0*'"“'^  sin[(N.i)iy] 

S-cos9,y^l^e  e  ^  ^  si„(m>) 

...  „  nr  I’/  *»»/ e-l®  sin(™“)  sin[(N.i)lfi] 

I2n  VR  0^  sinm 


(4-12) 
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The  average  power  pattern  for  the  grating  with  no  errors  is  given  by: 


1  k(Wcos0if 

sin(^) 

2 

sin 

(N+O  — 

N  +  1  27cR 

sin  1 

kTu^ 

1  2  i 

1 

^2 


1  k(Wcos0i)^ 

■sin(^)" 

2 

r  .  r,  vkTu  1 
sm  (N+O— 

N  +  1  27tR 

.  (^) . 

sin  (tS) 

(4-13) 


where 


Io  = 


k(Wcos0i)^ 

27iR 


(4-14) 


is  the  power  intensity  for  a  single  slit  or  strip  with  u=0.  With  the  exception  of  the  1/(N+1) 
factor.  Equation  (4-13)  is  the  familiar  Fraunhofer  diffraction  pattern  of  an  array  of  N+1 
slits  found  in  optics  textbooks  (10:410). 

To  find  the  average  power  pattern  for  the  infinite  strip  grating,  use  Equation  (2-7): 


Grating  with  Errors 

Now  let  the  width  of  the  strips  and  slits  be  defined  as  in  Equation  (3-53).  The  far 
zone  scattered  fields  in  Equation  (4-10)  become: 
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(4-16) 


=cosei-|^e'/^  f" 

llK  VR 

_  jCOSGj  £_'_'V  (giku^n  _  gjkuXn  \ 

uV27ik  Vr  „ 


Again  let  n  range  from  0  to  N.  The  average  power  pattern  is  given  by: 


S^(u)  = 


1  COS^  0; 


N  + 1  27tku^R 


t  n=0  m=0 


Let 


N  N 


d=ES(  gjkuxn  _  gjkuxn  ^^g-jkut„  _  g-jkux„  ^ 

n=0  m=0 
N  N 

_  ^  ^  j’gjku(t„-x„,)  _^^jku(x„-Xn,)  _  gjku(t„-x„)  _  gjku(x„-x„) 
n=0 111=0 


(4-18) 


Separating  the  terms  for  which  n=m,  this  becomes: 


D  =  (N  +  l)+(N  +  l)-^e^'‘"^'" 


n=0 


n=0 


f  N  N-n 


(4-19) 


+2Re  ELh 

[n=l  m=0 


jku(xj|+jii  XjTi,)  ^  gj^^(^n+ni  ^iij)  _ ^j^^(^n+m  ^m)  _ pi^^(^n+in  ^m) 


or 
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(4-20) 


D  =  2(N-H)-2Re 


/"  N  ^ 

^gMx„-xJ 


Vn=0 


-i-2Re< 


N  N-n  - 

IK' 

.  n=l  m=0 


jku(xn+in-Xm)  ^gjku(x„+„-x„)  _ ^ jku(x„+„-Xm)  _  p,jku(x„+„-x„ 


Now  the  following  relations  are  true: 


'Cn-Xn=W„ 


(4-21) 


'Cn+m-'^m  =am+Wm+i  +am+i  +W^+2+-”+an+m-l  +W„+m 
n+m-1 


(4-22) 


‘n+m  ^m=^m+  +'^m+l+  ^m+l  +’ ' '  +W„+m-l  +  a„+m-l 


n+m-1 


k=m 


(4-23) 


'tn+m  -Xm  =  +^m+—+'^n+m-l  +an+m-l  +'^n+m 

n+m~l 


=  Wn+m+  X^^k+wJ 


k=m 


(4-24) 


^n+m  ”’"^m  “•"^m+l  *^^m+l“*  n+m-1  "^^n+m-l 

n+m-1 

=  am+  S(^k+wJ 

k=m+l 


(4-25) 


Also  define  the  following  characteristic  functions: 


©(u)  =  E{ej'‘“"} 
0(u)  =  E{e''‘“"} 


(4-26) 


+6 


For  the  types  of  errors  given  by  Equation  (3-53),  these  become: 


0(u)  =  E{e»"}  =  =  ef‘"’*'E{e«} 

®(u)  =  E{e'‘"}  =  =e'“E{e''“} 

Then 


f  N 


N  N-n 


E{D}  =  2(N  +  l)-2Re  ^0  -i-2Rej  [a>"0" +0"©" -O"©"""' 


Vn=0 


^  n=l  ni=0 


=  2(N  +  l)-2(N  +  l)Re{0] 


f  N 

+  2Re]^(N  +  l-n)[o"©"  -hO"0”  -<D”©"^' 
ln=l 


(4-28) 


or 


E{d}  =  2(N-H) 


1  -  Re{©}  +  Rei  ^(^>©)"  [2  -  ©  -  )^] 


Ln=l 


-2Rel|;n(4>0)"[2-e-}4 

U=i 


(4-29) 


Use  Equations  (2-18)  and  (2-19)  to  obtain: 


E{d}  =  2(N-H) 


1  -  Re{©}  -I-  Re 


(<!>©)_ 

i-(a>©) 


[2-@->^ 


-2  Re 


r(a>©)-»-N(o©)^''^  -(n-h)(o©)’^^' 

[2-0- 1/1 . 

[  [i-(a)©)]^ 

r  /©J 

(4-30) 


Substituting  Equation  (4-30)  into  Equation  (4-17),  the  average  power  pattern  becomes: 
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^  COS^  0; 


„f/  COS 

1^  <“^l  =1^ 


1-Re{0}  +  Re 


(a>0)-(o0)^"‘ 

1-(O0) 


2-0-  I 


N  +  1 


Re-^ 


(00)  +  N(O0)  -  (N  +  l)(O0) 


\N+1 


(4-31) 


or  using  Equation  (4-14), 


S'(u)|  =Io 


(kuW) 


l-Re{0}  +  ReJ 


(0)0) -(O0f^' 

1-(O)0) 


[2-0- 


='0j 


N  +  1 


■Re<^ 


(<D0)+N(O0)^^^  -(N  +  1)(O0) 

[l-((I>0)f 


(4-32) 


Note  the  similarity  of  Equation  (4-32)  to  Equation  (2-20).  Substituting  u— 0  into  the  first 
equation  in  Equation  (4-16)  yields: 


N  +  1  27rR 


w„w 


* 

m 


^:^[(N  +  1)e{w'}+(n"  +n)eMw}] 

i[E{w’}  +  NEHw}] 


(4-33) 


or 


|sf(o)|^  =:^[(N  +  1)W^  +  e{4^}]  (4-34) 
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For  the  case  of  an  infinite  number  of  slits  or  strips,  i.e.  N  — >  oo.  Equations  (2-22)  and 
(2-23)  can  be  utilized  to  obtain; 


=  I 


"  (kuW)^ 


l-Re{0}  +  ReJ 


(<D0) 

1-(O0) 


2-0-  I 


0j 


(4-35) 


To  examine  the  behavior  at  u  =  0,  again  let  u  ^  }x  +  joo  and  use  the  first  order  terms  of 
the  following  Taylor  series  expansions; 


0(p,  +  j(i))»  [p3(a)[l-fflka  +  jk;|ia]da 
=  l-(oka  + jkpa 

0(ti  +  j©)«  [p,^(w)[l-®kw  +  jkpw]dw 
Jo 

=  1  -  ©kw  +  jk|j.w 

With  the  types  of  errors  given  by  Equations  (3-53)  and  (3-54), 

a  =  A  +  E{8}  =  A 
w  =  W  +  E{^}  =  W 


(4-36) 


(4-37) 


and  Equation  (4-36)  becomes; 

<l)((j.  +  j©  )  «  1  -  ©kA  +  jkpA  « 

0(|a  +  j© ) «  I  -  ©kW  +  jkpW  « 

for  small  and  ©.  Substitution  into  Equation  (4-35)  yields; 


(4-38) 
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S^(|a  +  jo) 


=  1. 


k(|^  +  j®)' 


W 


l2 


+Re<^ 


0  pjk(^l+j«>)w  -jk(n+ja))w'| 

► 

j  _  ^gjk(^+j®)Agjk(n+ja))wJ 

J 

(4-39) 


or 


S‘'(|i  +  j(o)|  =Io 


k(^i  +  j©fw]' 


-[l-cos[k(|x  + jco)w] 


+Re( 


^gjk(^+j«>)(A+W)j 
jk(n+jo))(A+W) 


1-e 


[2  -  2  cos[k((i  +  jo)  )w]] 


(4-40) 


which  becomes: 


S‘‘(^  +  j©)|  =Io 


-2<  sin 


k(fi  + j0)W 


2 


[k(^  +  j(ofw]  I  L  ^ 

^gjk(M+j“)(A+W)j 


1  +  Re( 


1-e 


jk(^+j(B)(A+W) 


(4-41) 


or 


S‘'(n  +  jco)|  =Io 


fsm[ 


k(n-hj<o)W 


k(n+j(»)W 


1- 


jk(n+j(o)T 


1-e 


jk(n+jio)T 


(4-42) 


Using  Equation  (2-3 1),  this  becomes  as  ®  0 ; 
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which  is  the  same  as  for  the  case  of  the  grating  with  no  errors. 


Grating  with  Uniformly  Distributed  Errors 

Suppose  8  is  uniformly  distributed  over  the  range;  -cA  <  e  <  cA  where  c  is  the 
maximum  spacing  error  expressed  as  a  percentage  of  the  desired  spacing  A.  The  PDF  of  e 
is  given  by: 


P.(e)^ 


1 


2cA 

0 


-  cA  <  8  <  cA 
otherwise 


The  characteristic  function  <I>  is  given  by; 


<D(u)  =  ej‘‘“^E{ej'‘‘“} 


eJkuA 

2cA 


^ikuA  sin(kucA) 
kucA 


(4-44) 


(4-45) 


Now  let  ^  be  uniformly  distributed  over:  -bW  <  ^  <  bW  where  b  is  the  maximum  width 
error  expressed  as  a  percentage  of  the  desired  width  W.  The  PDF  of  ^  is  given  by: 


1 

2bw 


[0 


-  bw  <  ^  <  bw 
otherwise 


(4-46) 


The  characteristic  function  0  is  given  by: 


(4-47) 


gjkuW 

2bW 


ikuw  sin(kubW) 
®  kubW 


Graphs  comparing  the  average  power  patterns,  normalized  by  dividing  by  Iq 
(denoted  by  I/Iq),  for  exact  gratings  and  gratings  with  errors  are  shown  below.  To  allow 
a  clearer  picture  of  the  effects  of  the  random  errors,  a  finite  number  of  slits/  strips  is  used 
for  most  of  the  graphs.  Plots  for  an  infinite  number  of  slits/  strips  are  included  to  show  the 
degradation  to  the  impulse  functions  of  the  ideal  grating.  The  following  notation  is  used 
in  the  graphs  to  indicate  the  values  of  the  appropriate  variables: 


u:  sin  Gj  -  sin  0 

N+1 :  Number  of  Slits  or  Strips. 

W/X:  Desired  Width  (Ratio  of  desired  slit  (TEz)  or  strip  (TMz)  width  (W)  to 
Wavelength) 

A/X:  Desired  Spacing  (Ratio  of  desired  strip  (TEz)  or  slit  (TMz)  width  (A)  to 
Wavelength) 

Uniform  Dist;  Values  for  grating  with  errors  computed  using  uniform  distributions  for 
the  errors 

%W  Tol:  The  maximum  value  allowed  of  the  error  in  width  (^)  expressed  as  a 

percentage  of  width  W.  This  corresponds  to  the  variable  "b"  in  Equation  (4-46) 

%A  Tol:  The  maximum  value  allowed  of  the  error  in  width  (e)  expressed  as  a  percentage 
of  width  A.  This  corresponds  to  the  variable  "c"  in  Equation  (4-44) 
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Figure  4.2  shows  the  average  power  pattern  for  an  exact  grating  (i.e.  a  grating 
with  no  errors)  with  a  finite  number  of  slits  or  strips.  The  plots  consist  of  a  series  of  large 
spikes  located  arcsin  (k/T)  apart  in  angle.  The  plot  clearly  shows  the  sine  envelope  of  the 
average  power  pattern  which  is  governed  by  the  w/T  ratio.  On  an  expanded  scale,  Figure 
4.3  shows  the  effects  of  increasing  the  total  number  of  slits  or  strips.  As  the  number  of 
shts/  strips  increases,  the  amplitudes  of  the  spikes  increase  and  while  the  widths  of  the 
spikes  decrease  and  the  pattern  begins  to  resemble  a  series  of  impulse  functions  as  given  in 
Equation  (4-15).  Figure  4.4  presents  the  same  pattern  with  the  values  normalized  so  that 
the  maximum  value  (at  u  =  0)  equals  one,  again  showing  the  sine  envelope  to  the  ratios  of 
the  main  spikes. 


Figure  4.2  Far  Field  Born  Approximation  -  Grating  with  no  Errors 
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re  4.3  Far  Field  Born  Approximation  -  Effects  of  Total  Number  of  Slits/ 
Strips  on  Grating  with  no  Errors 


Figure  4.4  Normalized  Far  Field  Born  Approximation  •  Effects  of  Total 
Number  of  Slits/  Strips  on  Grating  with  no  Errors 
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Figures  4.5  through  4.8  show  the  effects  of  errors  in  spacing  alone,  width  alone, 
and  both  spacing  and  width  respectively  for  a  finite  grating  with  desired  width  equal  to 
desired  spacing,  i.e.  W  =  A.  The  plots  show  the  overall  effect  is  to  lower  the  amplitude 
of  the  main  peaks  and  broaden  the  width  of  the  spikes.  The  values  between  the  main 
spikes  also  increase  slightly.  The  effects  become  more  pronounced  as  the  parameter  u  is 
increased.  A  comparison  of  Figure  4.5  to  Figure  4.6  shows  that  the  effects  of  error  in 
width  alone  are  on  the  order  of  the  errors  in  spacing  alone  for  a  given  tolerance  level. 
Figures  4.7  and  4.8  show  the  combined  effects  of  errors  in  both  spacing  and  width. 

Figures  4.9  through  4.12  show  the  effects  of  errors  in  spacing  alone,  width  alone, 
and  both  spacing  and  width  respectively  for  a  finite  grating  with  desired  width  less  than 
the  desired  spacing,  i.e.  W  <  A.  The  plots  again  show  the  overall  effect  is  to  lower  the 
amplitude  of  the  main  peaks  and  broaden  the  width  of  the  spikes.  A  comparison  of  Figure 
4.9  to  Figure  4.10  shows  that,  for  this  case,  the  effects  of  error  in  spacing  alone  are 
greater  than  the  effects  due  to  errors  in  width  alone  for  a  given  percent  tolerance  level.  In 
terms  of  absolute  tolerance  in  wavelengths,  however,  the  errors  are  similar,  i!e.  the  %A 
Tol  =  5  errors  in  Figure  4.9  are  similar  to  the  %W  Tol  =15  errors  in  Figure  4.10.  Figures 
4. 1 1  and  4.12  show  the  combined  effects  of  errors  in  both  spacing  and  width. 

Figures  4.13  through  Figure  4.16  show  the  effects  of  errors  in  spacing  alone, 
width  alone,  and  both  spacing  and  width  respectively  for  a  finite  grating  with  desired 
width  greater  than  the  desired  spacing,  i.e.  W  >  A.  A  comparison  of  Figure  4.13  to  Figure 
4.14  shows  that,  opposite  to  the  previous  case,  the  effects  of  error  in  width  alone  are 
greater  than  the  effects  due  to  errors  in  spacing  alone  for  a  given  percent  tolerance  level. 
In  terms  of  absolute  tolerance  in  wavelengths,  however,  the  errors  are  similar,  i.e.  the  %A 
Tol  =15  errors  in  Figure  4.13  are  similar  to  the  %W  Tol  =  5  errors  in  Figure  4.14. 
Figures  4. 15  and  4.16  show  the  combined  effects  of  errors  in  both  spacing  and  width. 
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gure  4.5  Far  Field  Born  Approximation  -  Effects  of  Spacing  Errors  on 
Grating  with  Desired  Width  Equal  to  Desired  Spacing 


Figure  4.6  Far  Field  Born  Approximation  -  Effects  of  Width  Errors  on 
Grating  with  Desired  Width  Equal  to  Desired  Spacing 
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Average  Power  (I/Io) 


Figure  4.7  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Equal  to  Desired  Spacing 


Figure  4.8  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Equal  to  Desired  Spacing 
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igure  4.9  Far  Field  Born  Approximation  -  Effects  of  Spacing  Errors  on 
Grating  with  Desired  Width  Less  Than  Desired  Spacing 
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Figure  4.10  Far  Field  Born  Approximation  -  Effects  of  Width  Errors  on 
Grating  with  Desired  Width  Less  Than  Desired  Spacing 
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;ure  4.11  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Less  Than  Desired  Spacing 
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Figure  4.12  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Less  Than  Desired  Spacing 
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Figure  4.13  Far  Field  Born  Approximation  -  Effects  of  Spacing  Errors  on 
Grating  with  Desired  Width  Greater  Than  Desired  Spacing 
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Figure  4.14  Far  Field  Born  Approximation  -  Effects  of  Width  Errors  on 
Grating  with  Desired  Width  Greater  Than  Desired  Spacing 
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gure  4.15  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Greater  Than  Desired  Spacing 


Figure  4.16  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Desired  Width  Greater  Than  Desired  Spacing 
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Figure  4.17  shows  the  average  power  pattern  for  a  grating  with  an  infinite  number 
of  slits/  strips  for  varying  amounts  of  equal  width  and  spacing  tolerances.  The  effects  due 
to  increasing  the  amount  of  errors  follow  the  same  trends  as  for  the  finite  cases.  Although 
not  shown  on  the  graph,  the  patterns  still  contain  an  impulse  at  the  origin  (u  =  0)  which  is 
not  plotted  due  to  its  infinite  amplitude.  Using  a  ray  analysis  approach,  the  presence  of 
the  impulse  can  be  regarded  as  being  due  to  the  fact  that  all  rays  travel  equal  optical  path 
lengths  for  u  =  0.  The  graph  also  shows  significant  reduction  in  amplitude  and  an 
increase  in  beam  widths  as  compared  to  a  grating  with  no  errors.  The  figure  also  shows 
non-zero  values  at  locations  between  the  main  spikes.  These  effects  increase  as  the 
parameter  u  is  increased,  essentially  washing  out  the  grating  lobe  pattern  of  the  grating 
with  no  errors. 

Figure  4.18  shows  the  progression  of  increasing  the  number  of  slits  or  strips  from  a 
small  finite  value  up  to  an  infinite  number.  Here  the  build  up  to  the  impulse  at  the  origin  is 
seen  as  well  as  the  to  the  finite  values  at  the  other  main  spike  locations. 

For  completeness  and  later  reference.  Figures  4.19  and  4.20  show  the  average 
power  pattern  for  low  frequency  scattering,  i.e.  T/A.  <1.  Figure  4.19  presents  the  results 
for  a  grating  with  a  finite  number  of  slits  or  strips.  The  main  effect  of  the  errors  is  to  raise 
the  value  of  the  nulls  in  the  floor  of  the  pattern  while  the  main  lobe  is  not  significantly 
affected.  Figure  4.20  presents  the  results  for  a  grating  with  an  infinite  number  of  slits  or 
strips. 
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ure4.17  Far  Field  Born  Approximation  -  Effects  of  Width  and  Spacing 
Errors  on  Grating  with  Infinite  Number  of  Slits  or  Strips 
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Figure  4.18  Far  Field  Born  Approximation  -  Effects  of  Total  Number  of 
Slits/  Strips  on  Grating  with  Errors 
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Figure  4.19  Far  Field  Born  Approximation  -  Effects  of  Errors  on  Grating 
with  Finite  Number  of  Slits  or  Strips,  Low  Frequency  Case 


Figure  4.20  Far  Field  Born  Approximation  -  Effects  of  Errors  on  Grating 
with  Infinite  Number  of  Slits  or  Strips,  Low  Frequency  Case 
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Grating  with  Cosine  Function  Distributed  Errors 

The  uniformly  distributed  errors  can  be  regarded  as  representing  a  worst  case 
distribution  of  the  errors  and  may  not  be  representative  of  typical  manufacturing  errors.  A 
PDF  based  upon  a  cosine  function  is  more  representative  of  actual  errors.  Let  the  PDF  of 
£  be  given  by: 


P,(e)  =  <{2cAL 
0 


1 


1  +  cosi 


-cA<e<cA 

otherwise 


(4-48) 


This  PDF  distributes  most  of  the  errors  around  zero  while  the  density  tapers  off  to  zero  at 
e  =  ±cA .  The  characteristic  function  O  is  given  by: 


0(u)  =  ej‘‘“^E{ej*‘“"} 

.cA 


2cA  •'-cA 


r 


jkue 


1  +  cosl 


(S)]de 


(4-49) 


Using  the  integral : 


J  e“  cosbx  dx 


e“  (a  cos  bx  -H  b  sin  bx) 


(4-50) 


the  characteristic  function  becomes: 


sin(kucA) 

kucA 


(4-51) 


Use  of  L'Hopital's  rule  yields: 
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n 


(4-52) 


0(U)  = 


for  u  = 

kcA 


Similarly,  let  the  PDF  of  ^  be  given  by: 


1 

2bW 

0 


f 

1  +  cos 

V 


-bW<|<bW 

otherwise 


(4-53) 


The  characteristic  function  @  is  given  by: 

©(u)  =  ej*'“'^E{eJ’''^} 

sin(kubW) 
j_^kubw^2  kubW 

and 

©(u)  =  ^e^  for  u  =  -^ 
2  kbW 


(4-54) 


(4-55) 


Graphs  showing  the  average  power  pattern  of  gratings  with  errors  using  the  cosine  PDF 
are  shown  below.  Figures  4.21  and  4.22  show  the  pattern  for  a  grating  with  a  finite 
number  of  sUts/  strips  while  Figure  4.23  shows  the  pattern  for  a  grating  with  an  infinite 
number  of  slits/  strips.  The  overall  effects  are  similar  to  the  gratings  with  the  uniformly 
distributed  PDFs.  As  expected,  for  a  given  maximum  error  (i.e.  %  Tol),  the  errors  in 
spacing  and  width  have  less  of  an  effect  upon  the  average  power  pattern  of  a  grating  with 
cosine  distributed  errors  than  to  a  grating  with  uniformly  distributed  errors. 
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Figure  4.21  Far  Field  Born  Approximation  -  Cosine  PDF,  Finite  Number  of 

Slits/Strips 


Figure  4.22  Far  Field  Born  Approximation  -  Cosine  PDF,  Finite  Number  of 

Slits/Strips 
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Figure  4.23  Far  Field  Born  Approximation  -  Cosine  PDF,  Infinite  Number  of 

Slits/Strips 
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V.  Edge  Condition  Approach 


Overview 

This  chapter  presents  results  obtained  when  the  unknown  current  in  the  integral 
field  equation  is  replaced  by  an  approximation  which  satisfies  the  edge  condition,  i.e,  the 
current  behaves  as  p”^'^  where  p  is  the  distance  from  the  edge  (11:577-579).  The 
approximation  is  first  obtained  for  the  case  of  a  grating  with  no  errors  and  an  infinite 
number  of  slits/strips  and  then  modified  for  the  case  of  a  grating  with  width  and  spacing 
errors.  This  approach  is  especially  valid  in  the  low  frequency  region  where  the  wavelength 
is  much  larger  than  the  widths  of  all  the  slits  and  strips. 


Far  Field  Pattern  for  Grating  with  No  Errors 

Setting  the  origin  at  the  mid  point  of  the  first  slit  or  strip,  i.e.  1/2  (Xq  +  Tq)  =  0, 
the  scattered  fields  for  a  grating  with  no  errors  are  given  by  (Equation  (3-44)); 

S  =  Is(x')H®(k7(x-x')^y')dx'  (5-1) 

2  JnT-W/2  * 

n 

which  becomes  in  the  far  field; 


S'  = 


27t 


e"-''*  v-'  (•nT+W/2 


(5-2) 


The  unknown  current  Ig  can  be  approximated  by  the  following  expression  (1 1 : 192): 


Isn(x')  = 


(5-3) 


5-1 


where  the  subscript  "n"  refers  to  the  current  in  the  n^h  slit  or  strip  and  "C"  is  a  constant. 
The  value  for  the  constant  "C"  can  be  obtained  using  a  modal  expansion  for  the  scattered 
field  as  provided  in  Appendix  B.  This  results  in; 


C  = 


1 


P.T 


(5-4) 


where 


k2-(ksin0i+^f 


iK 


(5-5) 


Figure  5.1  compares  the  reflection  and  transmission  coefficients  as  a  function  of  the  TlX 
ratio  computed  using  this  approximation  to  the  coefficients  obtained  from  an  exact 
solution  using  the  Wiener-Hopf  method  (29:267-281).  The  coefficients  are  calculated  for 
a  grating  with  no  errors,  an  infinite  number  of  slits/strips,  normal  incidence  TEz  mode,  and 
the  width  equal  to  the  spacing,  i.e.  W  =  A  =  T/2.  For  the  range  of  T/X  ratios  shown, 
corresponding  to  only  one  propagating  scattered  mode,  the  approximation  agrees  very 
well  with  the  exact  solution. 


Substituting  Equation  (5-3)  into  Equation  (5-2)  results  in: 


sf  = 


k  •’%  e  Y>  j-nT+W/S 
2^  ® 


Ce 


jkx'sinBj 


[('%)^-(x'-nT)2] 


V2 


g.jkx'sinej^, 


(5-6) 


Using  the  change  of  variables: 


x”  =  x'-nT 


(5-7) 
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Figure  5.1  Comparison  of  Approximation  Using  the  Edge  Condition  to 
Exact  Solution  Obtained  via  the  Wiener-Hopf  Method 


Equation  (5-6)  becomes: 


*  ^  dx 


271 


Vr  4-J- 


W/2 


^  fk  ^y  e 
=  C-  —  e  — p=- 
)  27c  VR 


Jkux" 


(5-8) 


— _ 

.  J-w/= 


dx' 


The  integral  is  simplified  in  Appendix  C  with  the  following  results: 


fW/2 

J-W/2 


jkux" 


[(W^)2_(x")2] 


H 


dx"  =  7tJo(^) 


(5-9) 
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Which  leads  to: 


\cTr  j^/  e 

Sf  =  C,He  ^  V- 

V  2  ^/R 


(5-10) 


The  average  power  pattern  for  the  grating  with  no  errors  is  given  by; 


sL»(uf=^^[j«(^)r 


‘  .  kTu 


(5-11) 


where 


C^kTr 


(5-12) 


To  find  the  average  power  pattern  for  the  infinite  strip  grating,  use  Equation  (2-7): 


S[^(u)f=5^I[j.(^)r8(u 


(5-13) 


Grating  with  Errors 

Now  let  the  width  of  the  strips  and  slits  be  defined  as  in  Equation  (3-53).  Using 


Xn+^n+W„ 


Xn+'Cn-Wn 


(5-14) 


5-4 


The  far  zone  scattered  field  becomes: 


S' 


e-jkR 

Vr 


SI 


Xn+Tn+Wn 
2 


Xn-«-^n-Wn 
2 


I,(x')e-j'“'*‘"®dx' 


(5-15) 


The  approximation  for  the  unknown  scattering  current  given  in  Equation  (5-3)  is  modified 
to  retain  the  appropriate  edge  condition  behavior  as  follows: 


Isn(x')  =  r 


Ce 


jkx'sinGj 


H 


(5-16) 


Substituting  this  into  Equation  (5-15),  the  scattered  field  becomes: 


FT  in/  «-jkR  Xn-nn->wn  T  ~  (  \\2T'^ 


Using  the  change  of  variables: 


Y  4-  T 


(5-18) 


The  scattered  field  becomes: 


-(x"f 


dx” 


(5-19) 
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Using  Equation  (5-9),  this  becomes: 


S^=C 


e-jkR 

Vr 


(5-20) 


Letting  n  range  from  0  to  N,  the  average  power  pattern  is  given  by: 


C^Ttk 

(N-fl)2R 


N  N 


Eli 

ln=0  m=0 


j¥(xn+^n)T  (kuwn)  -jyUm+^m)^ 

e  Jol-2^je  ^  Jo 


(5-21) 


Separating  terms  for  which  n=m, 


C^Tik 
(N  +  1)2R 


n=l  m=0 


(5-22) 


or 


C^Ttk 


2R 


[E{[j„(l=)f} 


+  -^Re| 
N  +  1 


N  N-n 

E|II« 

^  n=l  m=0 


+T 

J  2  *’n+m  ’’m 


(5-23) 


Let 


N  N-n 


B  =  IIe 


j  2  (^n+m'*"^n+m 


n=l  m=0 


(5-24) 
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From  Equations  (4-22)  and  (4-23),  the  following  is  true: 


Xn+m-Xm+'t„+„,-X:„=  X(^k+Wk)+ 

k=m  k=m 

n+m-1 

=  2a„+w„+w„+„-H2  ^(ak+wj 


(5-25) 


Using  the  characteristic  functions: 


0(u)  =  Ejej'^'*'" }  =  =  ej‘'“'^E{ej''“^} 

0(u)  =  E{ej''“"}  =  E{ej''“^^^"^}  =  ej‘'“^E{ej‘'“} 


(5-26) 


Then 


E{b}  =  X  Jo  W 


n=l  m=0 


E^  I  e^  Jo(^)  [X  J  “  n)(©a>)"  ©“^ 


[  kuw  I 

=  E^\ e^ Jo{-‘t)  “ S n(©0)" 


(5-27) 


Use  Equations  (2-18)  and  (2-19)  to  obtain: 


+n(©o)^^^^^  -  (n + l)(©o/^'"^ 

[i-©a)f 


(5-28) 


5-7 


Inserting  Equation  (5-28)  into  Equation  (5-23),  the  average  power  pattern  becomes: 


S'(u)f  =-^[E{[jo(^)f  }  +  2Re  E2p'rj^(k^)|0-i 


-2Re 


;kuw 


1-0(D 

+N(0O)^^‘'^^  -  (N  + 


(N  +  1)[1-0O] 


.n  (5-29) 


Substituting  u=0  into  Equation  (5-21)  yields: 

|s'(0)f  =-^(N  +  1)  =  Io(N  +  1)  (5-30) 

'  2R 


For  the  case  of  an  infinite  number  of  slits  or  strips,  i.e.  N  — >  oo ,  Equations  (2-22)  and 
(2-23)  can  be  utilized  to  obtain: 


(5-31) 


To  examine  the  behavior  at  u  =  0,  again  let  u  ^  |i  +  j©  and  use  the  first  order  terms  of 
the  following  Taylor  series  expansions  for  u  +  j®  «  0: 


T  ^  I 

°L  2  J~ 


(5-32) 


.k(^+j®) 

e^  2 


w 


«1- 


akw 

~2~ 


(5-33) 


which  leads  to: 


5-8 


(5-34) 


.k(u+j(»)w 


k(^+jco)w 

2 


cokw 


and 


.k(fA+jCD)w 


Jr 


k(^+j(o)w 


Pw(w) 


cokw  .  kuw 


dw 


cokw 

2 


(5-35) 


(5-36) 


With  the  types  of  errors  given  by  Equations  (3-53)  and  (3-54), 


a  =  A  +  e{s}  =  A 
w  =  W  +  e{^}  =  W 


(5-37) 


Equation  (5-35)  becomes: 


.k(n-l-jo))w 
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for  small  \i  and  ©.  Using  Equation  (4-38), 


0(p  -f  j®  )(l)(p  +  j© ) « 


(5-39) 


Substitution  into  Equation  (5-31)  yields: 
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or 
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Using  Equation  (2-31),  as  ©  — >  0  this  becomes: 


S‘'(|a  +  jo)|  = 


^2^5(h)  =  ^6(h) 
2R  T  ^  T 


for  =  0 


(5-42) 


Grating  with  Uniformly  Distributed  Errors 

Let  8  and  ^  have  uniform  PDFs  as  given  in  Equation  (4-44)  and  Equation  (4-46). 
The  corresponding  characteristic  functions  O  and  0  are  given  by  Equation  (4-45)  and 
Equation  (4-47).  Graphs  comparing  the  average  power  patterns  (normalized  by  dividing 
by  Iq)  for  exact  gratings  and  gratings  with  errors  are  shown  below.  The  majority  of  the 
graphs  present  results  for  low  frequency  scattering  i.e.  ?)>T  since  the  edge  condition 
approach  is  more  accurate  in  this  regime.  The  expected  values  for  terms  involving  the 
Bessel  function  (Jq)  were  computed  numerically.  To  allow  comparison  of  gratings  with 
no  errors  to  gratings  with  errors,  a  number  of  graphs  were  computed  using  a  finite  number 
of  slits  or  strips.  This  permits  the  graphs  to  show  the  progressive  changes  in  the  average 
power  pattern  of  a  grating  with  errors  to  a  grating  with  no  errors. 
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Figure  5.2  shows  the  effects  of  increasing  the  total  number  of  slits  or  strips  for  a 
grating  with  no  errors  and  T  =  0. 8A, .  As  expected,  as  the  number  of  slits  or  strips 
increases,  the  pattern  begins  to  resemble  an  impulse  function  at  the  origin.  This  is 
consistent  with  only  one  propagating  mode  being  scattered  from  a  grating  with  no  errors 
and  an  infinite  number  of  slits  or  strips  when  A,  >  T. 


u 


Figure  5.2  Far  Field  Edge  Condition  Approximation  -  Effects  of  Total 
Number  of  Slits/  Strips  on  Grating  with  no  Errors 

Figures  5.3  to  5.11  show  the  effects  of  errors  in  spacing  alone,  width  alone,  and 
both  spacing  and  width  for  finite  gratings  with  various  ratios  of  desired  width  to  desired 
spacing.  The  desired  period  for  the  graphs  is  T  =  0.8A.  The  plots  show  the  overall 
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effects  are  similar  to  the  effects  seen  using  the  zero  order  approximation.  Most 
significantly,  the  main  lobe  is  not  significantly  degraded  and  the  nulls  in  the  floor  are 
raised. 

Figures  5.3  to  5.5  show  the  effects  of  errors  in  spacing  alone,  width  alone,  and 
both  spacing  and  width  for  a  finite  grating  with  desired  width  equal  to  desired  spacing,  i.e. 
W  =  A.  A  comparison  of  Figure  5.3  to  Figure  5.4  shows  that,  as  in  the  results  obtained 
using  the  zero  order  approximation,  the  effects  of  error  in  width  alone  are  on  the  order  of 
the  errors  in  spacing  alone  for  a  given  tolerance  level.  Figure  5.5  shows  the  combined 
effects  of  errors  in  both  spacing  and  width.  Note  the  simularity  of  this  graph  to  Figure 
4.19,  obtained  using  the  zero  order  approach. 

Figures  5.6  through  5.8  show  the  effects  of  errors  in  spacing  alone,  width  alone, 
and  both  spacing  and  width  respectively  for  a  finite  grating  with  desired  width  less  than 
the  desired  spacing,  i.e.  W  <  A.  A  comparison  of  Figure  5.6  to  Figure  5.7  shows  that,  for 
this  case,  the  effects  of  error  in  spacing  alone  are  again  greater  than  the  effects  due  to 
errors  in  width  alone  for  a  given  percent  tolerance  level.  Figure  5.8  shows  the  combined 
effects  of  errors  in  both  spacing  and  width. 

Figures  5.9  through  5.11  show  the  effects  of  errors  in  spacing  alone,  width  alone, 
and  both  spacing  and  width  respectively  for  a  finite  grating  with  desired  width  greater  than 
the  desired  spacing,  i.e.  W  >  A.  A  comparison  of  Figure  5.9  to  Figure  5.10  shows  that, 
opposite  to  the  previous  case,  the  effects  of  error  in  width  alone  are  again  greater  than  the 
effects  due  to  errors  in  spacing  alone  for  a  given  percent  tolerance  level.  Figure  5.11 
shows  the  combined  effects  of  errors  in  both  spacing  and  width. 
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Average  Power  (I/Io) 


Figure  5.3  Far  Field  Edge  Condition  Approximation  -  Effects  of  Spacing 
Errors  on  Grating  with  Desired  Width  Equal  to  Desired  Spacing 


Figure  5.4  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width 
Errors  on  Grating  with  Desired  Width  Equal  to  Desired  Spacing 
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ure  5.5  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width  and 
>acing  Errors  on  Grating  with  Desired  Width  Equal  to  Desired  Spacing 
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Figure  5.6  Far  Field  Edge  Condition  Approximation  -  Effects  of  Spacing 
Errors  on  Grating  with  Desired  Width  Less  Than  Desired  Spacing 

5-14 


Figure  5.7  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width 
Errors  on  Grating  with  Desired  Width  Less  Than  Desired  Spacing 


Figure  5.8  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width  and 
Spacing  Errors  on  Grating  with  Desired  Width  Less  Than  Desired  Spacing 

5-15 


Figure  5.9  Far  Field  Edge  Condition  Approximation  -  Effects  of  Spacing 
Errors  on  Grating  with  Desired  Width  Greater  Than  Desired  Spacing 


Figure  5.10  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width 
Errors  on  Grating  with  Desired  Width  Greater  Than  Desired  Spacing 
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Figure  5.11  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width  and 
Spacing  Errors  on  Grating  with  Desired  Width  Greater  Than  Desired  Spacing 

Figure  5.12  shows  the  average  power  pattern  for  a  grating  with  an  infinite  number 
of  slits/  strips  for  varying  amounts  of  equal  width  and  spacing  tolerances  and  a  desired 
period  T  =  0.8X .  Although  not  shown  on  the  graph,  the  patterns  stiU  contain  an  impulse 
at  the  origin  (u  =  0)  which  is  not  plotted  due  to  its  infinite  amplitude.  Whereas  a  grating 
with  no  errors  would  have  a  power  pattern  equal  to  zero  for  u  0 ,  the  graph  shows  small 
but  non-zero  values  for  u  0 .  Note  the  simularity  of  this  graph  to  Figure  4.20. 

Figure  5.13  shows  the  average  power  pattern  for  a  grating  with  an  infinite  number 
of  slits/  strips  for  varying  amounts  of  equal  width  and  spacing  tolerances  and  a  desired 
period  T  =  20X, .  Although  not  shown  on  the  graph,  the  patterns  still  contain  an  impulse  at 
the  origin  (u  =  0)  which  is  not  plotted  due  to  its  infinite  amplitude.  The  graph  is  very 
similar  to  Figure  4.17  which  has  the  same  width  and  spacing  values. 
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re  5.12  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width  and 
Spacing  Errors  on  Grating  with  Infinite  Number  of  Slits  or  Strips 


1000 


0.0001  ■I — ' — ' — ' — ' — ^ — ' — ' — ' — '  '  '  '  '  '  '  '  '  ^ 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

u 

_ N+l=Infinite  W/l=\0  A/Ti^lO  Uniform  Dist 

Figure  5.13  Far  Field  Edge  Condition  Approximation  -  Effects  of  Width  and 
Spacing  Errors  on  Grating  with  Infinite  Number  of  Slits  or  Strips 
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Grating  with  Cosine  Function  Distributed  Errors 

Let  e  and  ^  have  cosine  function  PDFs  as  given  in  Equation  (4-48)  and  Equation 
(4-53).  The  corresponding  characteristic  functions  O  and  0  are  given  in  Equation  (4-51) 
and  Equation  (4-54).  Graphs  showing  the  average  power  pattern  of  gratings  with  errors 
using  the  cosine  PDF  are  shown  below.  Figure  5.14  shows  the  pattern  for  a  grating  with 
a  finite  number  of  slits/  strips  while  Figure  5.15  shows  the  pattern  for  a  grating  with  an 
infinite  number  of  slits/  strips.  The  overall  effects  are  similar  to  the  gratings  with  the 
uniformly  distributed  PDFs.  As  expected,  for  a  given  maximum  error  (i.e.  %  Tol),  the 
errors  in  spacing  and  width  have  less  of  an  effect  upon  the  average  power  pattern  of  a 
grating  with  cosine  distributed  errors  than  to  a  grating  with  uniformly  distributed  errors. 


Figure  5.14  Far  Field  Edge  Condition  Approximation  -  Cosine  PDF,  Finite 

Number  of  Slits/  Strips 


:e  Power  (I/Io) 


Figure  5.15  Far  Field  Edge  Condition  Approximation  -  Cosine  PDF,  Infinite 

Number  of  Slits/  Strips 
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VI.  Realizations  and  Validation 


Overview 

This  chapter  compares  the  statistical  average  power  pattern  computed  using  the 
expressions  developed  in  the  previous  chapters  to  the  average  of  a  number  of  trial 
realizations.  To  compute  the  average  power  pattern  for  a  realization,  the  errors  in  a 
grating  are  randomly  generated  using  the  uniform  or  cosine  function  PDFs.  For  all  cases, 
a  total  of  twenty  realizations  are  computed  for  each  grating.  The  average  power  pattern  is 
then  computed  for  each  grating  and  the  results  averaged  together.  Plots  showing  the 
standard  deviation  of  the  realizations  are  also  provided. 

Realizations  Using  the  Born  Approximation 

Figures  6.1  through  6.12  present  the  results  of  trial  realizations  using  the  Bom 
approximation.  Figure  6. 1  shows  a  typical  realization  of  the  average  power  pattern  for  a 
grating  with  N+l=  10  slits/  strips  and  uniformly  distributed  errors  with  maximum  width 
and  spacing  errors  equal  to  five  percent.  The  graph  also  shows  the  statistical  average 
power  pattern  computed  using  the  formulas  developed  in  Chapter  IV.  As  shown,  the 
realization  tends  to  follow  the  statistical  average.  Figure  6.2  shows  the  results  of 
averaging  the  power  pattern  for  twenty  trial  realizations.  Even  for  only  twenty 
realizations,  the  statistical  average  and  the  average  of  the  realizations  show  excellent 
agreement.  The  standard  deviation  and  the  standard  deviation/  mean  ratio  of  the  twenty 
realizations  are  shown  in  Figure  6.3.  The  standard  deviation  at  the  origin  (u  =  0)  is  small 
(-0.14).  For  values  of  u  greater  than  zero,  the  standard  deviation  resembles  the  average 
power  pattern.  In  fact,  the  ratio  of  the  standard  deviation  to  the  mean  shows  this  as  the 
curve  starts  out  at  a  small  value  at  u  =  0  and  then  oscillates  about  a  ratio  equal  to  one. 
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Figures  6.4  through  6.6  present  results  for  the  same  gratings  with  maximum  width 
and  spacing  errors  equal  to  ten  percent.  The  same  trends  are  apparent  and  the  statistical 
average  power  pattern  again  agrees  with  the  average  of  the  realizations.  The  number  of 
slits/  strips  is  increased  to  N+l=  100  for  Figures  6.7  through  6.9.  Again  the  statistical 
average  power  pattern  agrees  with  the  average  of  the  realizations  and  the  ratio  of  the 
standard  deviation  to  the  mean  starts  out  at  a  small  value  at  the  origin  and  then  oscillates 
about  a  ratio  equal  to  one.  This  behavior  of  the  standard  deviation  is  consistant  with  the 
expression  developed  in  Appendix  D  for  the  standard  deviation  of  a  grating  with  an 
infinite  number  of  slits/  strips. 

Finally,  Figures  6.10  through  6.12  present  the  results  for  a  grating  with  N+l=  10 
slits/  strips  with  errors  generated  using  the  cosine  function  PDF.  Again,  the  same  trends 
appear  as  for  the  gratings  with  the  uniformly  distributed  errors. 
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N+l=  10  W/X=  10  A/X-  10  %  W  Tol-  5  %  A  Tol=  5  Uniform  Dist 


Figure  6.1  Realization  of  Far  Field  Average  Power  Pattern  -  Born 
Approximation  with  Uniformly  Distributed  Errors 
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Std  Dev  or  Std  Dev/Mean 


Figure  6.2  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
Realizations  to  Statistical  Average  -  Born  Approximation  with  Uniformly 

Distributed  Errors 


Figure  6.3  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of 
Realizations  -  Born  Approximation  with  Uniformly  Distributed  Errors 


6-3 


Figure  6.4  Realization  of  Far  Field  Average  Power  Pattern  -  Born 
Approximation  with  Uniformly  Distributed  Errors 
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N+l=10  W/X-10  A/X.=  10  %WTol=10  %ATol=10  Uniform  Dist 


Figure  6.5  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
Realizations  to  Statistical  Average  -  Born  Approximation  with  Uniformly 

Distributed  Errors 
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Figure  6.6  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of 
tealizations  -  Born  Approximation  with  Uniformly  Distributed  Errors 
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Figure  6.7  Realization  of  Far  Field  Average  Power  Pattern  -  Born 
Approximation  with  Uniformly  Distributed  Errors 
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igure  6.8  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
ealizations  to  Statistical  Average  -  Born  Approximation  with  Uniformly 

Distributed  Errors 


Figure  6.9  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of 
Realizations  -  Born  Approximation  with  Uniformly  Distributed  Errors 
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N+l=  10  W/A.=  10  A/X=  10  %  W  Tol=  10  %  A  Tol=  10  Cosine  Dist 

Figure  6.10  Realization  of  Far  Field  Average  Power  Pattern  -  Born 
Approximation  with  Cosine  Distributed  Errors 
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Figure  6.11  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
Realizations  to  Statistical  Average  -  Born  Approximation  with  Cosine 

Distributed  Errors 
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Figure  6.12  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of 
Realizations  -  Born  Approximation  with  Cosine  Distributed  Errors 
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Realizations  Using  the  Edge  Condition  Approximation 

Figures  6.13  through  6.15  present  the  results  of  trial  realizations  using  the  Edge 
Condition  approximation.  Figure  6.13  shows  a  typical  realization  of  the  average  power 
pattern  for  a  grating  with  N+l=  10  slits/  strips  and  uniformly  distributed  errors  with 
maximum  width  and  spacing  errors  equal  to  ten  percent.  The  graph  also  shows  the 
statistical  average  power  pattern  computed  using  the  formulas  developed  in  Chapter  V. 
As  shown,  the  realization  tends  to  follow  the  statistical  average.  Figure  6.14  shows  the 
results  of  averaging  the  power  pattern  for  twenty  trial  realizations.  The  statistical  average 
and  the  average  of  the  realizations  show  excellent  agreement.  The  standard  deviation  and 
the  standard  deviation/  mean  ratio  are  shown  in  Figure  6.15.  As  a  consequence  of  using 
the  edge  condition  approximation,  the  standard  deviation  at  the  origin  is  zero.  For  values 
of  u  greater  than  zero,  the  standard  deviation  resembles  the  average  power  pattern.  The 
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ratio  of  the  standard  deviation  to  the  mean  follows  the  same  trend  has  for  the  previous 
cases  although  the  ratio  does  not  quite  reach  a  value  of  one  about  which  to  oscillate.  The 
number  of  slits/  strips  is  increased  to  N+l=  100  for  Figures  6.16  through  6.18.  Again  the 
statistical  average  power  pattern  agrees  with  the  average  of  the  realizations.  As  shown  in 
Figure  6.18,  the  ratio  of  the  standard  deviation  to  the  mean  does  reach  a  value  of  one 
about  which  to  oscillate.  Again,  the  behavior  of  the  standard  deviation  is  consistant  with 
the  expression  developed  in  Appendix  E  for  the  standard  deviation  of  a  grating  with  an 
infinite  number  of  slits/  strips. 


Figure  6.13  Realization  of  Far  Field  Power  Pattern  -  Edge  Condition 
Approximation  with  Uniformly  Distributed  Errors 
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Figure  6.14  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
Realizations  to  Statistical  Average  -  Edge  Condition  Approximation  with 

Uniformly  Distributed  Errors 


Figure  6.15  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of  Realizations  - 
Edge  Condition  Approximation  with  Uniformly  Distributed  Errors 
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Figure  6.16  Realization  of  Far  Field  Power  Pattern  -  Edge  Condition 
Approximation  with  Uniformly  Distributed  Errors 
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Figure  6.17  Far  Field  Average  Power  Pattern  -  Comparison  of  Average  of 
Realizations  to  Statistical  Average  -  Edge  Condition  Approximation  with 

Uniformly  Distributed  Errors 
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Figure  6.18  Far  Field  Average  Power  Pattern  -  Standard  Deviation  of 
Realizations  -  Edge  Condition  Approximation  with  Uniformly  Distributed 

Errors 
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VII.  Conclusions 


Overview 

This  chapter  summarizes  the  accomplishments  of  this  research  effort,  provides 
observations  and  conclusions  regarding  the  results,  and  suggests  several  areas  for  possible 
follow  on  investigation. 

Accomplishments 

The  first  accomplishment  of  this  effort  is  the  formulation  for  the  statistical  average 
power  pattern  of  a  linear  array  of  point  sources  with  errors  in  spacing.  This  is  equivalent 
to  the  average  power  pattern  of  a  linear  antenna  array.  Other  authors  have  investigated 
the  problem  of  the  random  array  by  treating  the  location  of  each  element  as  random.  The 
formulation  developed  in  this  work  treats  the  spacing  between  the  elements  as  random  or 
as  a  desired  spacing  plus  or  minus  a  random  error.  This  approach  is  more  practical  from  a 
manufacturing  standpoint  than  the  random  location  approach.  Although  the  development 
of  the  average  power  pattern  of  a  linear  array  of  point  sources  is  presented  to  provide  a 
basis  for  the  formulation  of  the  grating  problem,  its  applicability  to  antenna  array  analysis 
makes  it  a  noteworthy  accomplishment  in  its  own  right.  This  is  especially  true  when  one 
considers  the  ease  with  which  random  amplitude  and  phase  errors  could  be  added  to  the 
formulation. 

The  major  accomplishment  of  this  effort  is  the  development  of  two  formulations 
for  the  far-field  statistical  average  power  pattern  of  a  strip  grating  with  errors  in  the  widths 
of  the  strips  and  the  slits  as  compared  to  an  ideal,  exact  periodic  grating.  The  first 
formulation  utilizes  the  Bom  approximation  in  which  the  unknown  aperture  fields  are 
replaced  by  the  incident  fields.  The  second  formulation  utilizes  an  approximation  for  the 
unknown  fields  which  satisfies  the  edge  condition.  Approximations  for  the  scattered  fields 
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are  first  derived  using  PEC  surface  equivalence  for  a  TEz  polarized  plane  wave  incident 
upon  a  strip  grating  consisting  of  an  infinite  PEC  screen  cut  by  a  number  of  infinitely  long 
slits  (infinite  in  the  z-direction).  Babinet's  principle  is  then  used  to  obtain  approximations 
for  the  fields  obtained  when  a  TMz  polarized  plane  wave  is  incident  upon  the 
complementary  grating  formed  by  interchanging  the  slits  and  strips  of  the  original  grating, 
i.e.  a  number  of  infinitely  long  strips  in  free  space.  Expressions  for  the  average  power 
pattern  are  developed  in  terms  of  the  following  variables:  number  of  slits  or  strips,  desired 
width-  and  spacing-to-wavelength  ratios,  and  the  characteristic  functions  of  the  PDFs  of 
the  width  and  spacing  errors. 

The  formulation  of  the  far-field  statistical  average  power  pattern  is  unique  in  that  it 
is  based  upon  independent  but  identically  distributed  errors  in  the  widths  of  each  slit/  strip 
and  independent  but  identically  distributed  errors  in  the  spacings  between  them.  The  use 
of  the  independent  errors  leads  to  expressions  incorporating  geometric  series  involving  the 
characteristic  functions  of  the  width  and  spacing  PDFs  which  can  be  evaluated  in  closed 
form  even  for  an  infinite  number  of  terms.  This  allows  the  average  power  pattern  of  a 
grating  with  errors  to  be  calculated  for  gratings  with  an  infinite  as  well  as  a  finite  number 
of  strips  and  slits. 

Observations  and  Conclusions 

A  number  of  observations  and  conclusions  can  be  drawn  from  the  graphs  in 
Chapters  IV  and  V  which  show  the  far-field  average  power  patterns  of  gratings  with 
different  geometries.  Three  observations  are  immediately  apparent  when  looking  at  the 
patterns  of  gratings  with  a  finite  number  of  slits/  strips.  First  and  possibly  the  most 
significant,  the  main  lobe  of  a  grating  with  errors  is  not  significantly  different  from  that  of 
a  grating  with  no  errors.  Second,  the  grating  lobes  are  reduced  in  magnitude  and 
broadened  slightly.  Third,  the  level  of  the  sidelobes  between  the  grating  lobes  and  the 
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nulls  between  the  sidelobes  are  raised.  Furthermore,  the  effects  increase  as  the 
observation  angle  moves  away  from  the  main  lobe. 

The  graphs  also  show  the  relative  effects  of  the  errors  are  consistent  with  the  ratio 
of  the  desired  width  to  desired  spacing.  For  the  same  amount  of  maximum  percent  error, 
the  effects  of  spacing  errors  are  more  pronounced  on  a  grating  with  desired  spacing 
greater  than  desired  width  than  the  effects  of  width  errors.  Similarly,  for  the  same  amount 
of  maximum  percent  error,  the  effects  of  width  errors  are  more  pronounced  on  a  grating 
with  desired  width  greater  than  desired  spacing  than  the  effects  of  spacing  errors.  For  a 
grating  with  desired  width  equal  to  desired  spacing,  the  effects  of  errors  in  width  are  on 
the  order  of  the  effects  due  to  errors  in  spacing. 

The  same  observations  can  be  made  regarding  gratings  with  an  infinite  number  of 
slits/  strips.  For  a  periodic  grating  with  no  errors  and  an  infinite  number  of  slits  and  strips, 
the  power  pattern  consists  of  a  number  of  impulse  functions  located  arcsin(XyT)  apart  in 
angle.  These  impulse  functions  correspond  to  plane  waves.  As  shown  in  Chapters  IV  and 
V,  the  impulse  fimction  at  the  origin  (i.e.  the  main  lobe)  for  an  infinite  strip  grating  with 
errors  is  unchanged  from  that  of  a  grating  with  no  errors.  The  remaining  grating  lobes  are 
significantly  reduced  in  amplitude  and  broadened  in  width  as  evidenced  by  the  fact  they  are 
no  longer  impulse  fimctions.  This  means  that,  except  for  the  main  lobe,  the  scattered 
waves  are  no  longer  discrete  plane  waves.  The  power  patterns  of  the  infinite  strip  gratings 
also  show  non-zero  values  between  the  grating  lobes. 

It  is  interesting  to  note  that  the  average  power  pattern  for  a  grating  computed 
using  the  Bom  approximation  is  similar  to  the  average  power  pattern  for  a  grating 
computed  using  the  edge  condition  approximation.  This  leads  to  the  conclusion  that  the 
effects  would  be  similar  if  an  exact  expression  for  the  unknown  fields  could  be  obtained 
and  utilized. 
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The  power  patterns  of  a  number  of  strip  gratings  containing  random  errors  in 
spacing  and  width  were  computed  and  the  results  averaged  together.  As  shown  in 
Chapter  VI,  the  average  of  these  realizations  agrees  with  the  computed  statistical  average 
for  both  the  Bom  approximation  and  the  edge  condition  approximation.  This  result 
validates  the  formulations  developed  in  Chapters  IV  and  V  for  the  statistical  average 
power  pattern.  The  standard  deviation  of  the  realizations  was  also  computed  and  the 
results  agree  with  the  formulations  developed  in  Appendices  A,  D,  and  E  for  the  case  of  a 
grating  with  an  infinite  number  of  slits/  strips. 

The  behavior  of  the  standard  deviation  approaching  the  average  power  pattern  as 
the  number  of  slits/  strips  goes  to  infinity  is  quite  interesting.  As  a  possible  explanation, 
consider  the  case  of  the  array  of  point  sources.  By  invoking  the  central  limit  theorem  for 
probability,  the  PDF  of  the  array  factor  approaches  gaussian.  As  a  result,  the  PDF  of  the 
square  of  the  magnitude  of  the  array  factor,  i.e.  the  power  factor,  approaches  a  negative 
exponential.  One  property  of  a  process  with  a  negative  exponential  PDF  is  that  the 
standard  deviation  is  equal  to  the  mean. 

It  is  obvious  that  the  formulations  developed  in  this  work  can  be  used  to  determine 
the  average  power  pattern  for  a  strip  grating  for  a  given  distribution  of  errors  in  width  and 
spacing.  Conversely,  they  can  be  used  to  determine  allowed  tolerance  limits  to  meet 
specific  design  criteria.  Additionally,  the  formulation  developed  for  the  two-dimensional 
case  of  a  strip  grating  can  be  expanded  to  the  three-dimensional  case  of  square  or 
rectangular  apertures  in  a  PEC  plane  or  its  complement. 

Recommendations 

Several  areas  exist  where  follow  on  work  relating  to  this  research  effort  can  be 
investigated.  The  first  area  involves  an  extension  of  the  analysis  of  the  linear  array  of 
point  sources.  The  applicability  of  this  analysis  to  antenna  array  theory  makes  the  addition 
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of  random  phase  and  amplitude  errors  a  worthy  effort.  The  additional  effort  in  this  area 
could  be  strengthened  by  extending  the  analysis  to  two-dimensional  arrays. 

Second,  as  mentioned  above,  the  formulations  presented  for  the  average  power 
pattern  of  a  strip  grating  using  the  Born  and  edge  condition  approximations  can  be 
expanded  to  include  square  or  rectangular  apertures  on  a  PEC  plane  (or  its  complement) 
as  opposed  to  infinitely  long  slits  and  strips.  The  examinations  of  the  quadruple 
summations  found  in  the  Appendices  A,  D,  and  E  will  be  helpful  in  these  areas. 

A  third  area  of  possible  additional  effort  involves  the  investigation  into  a  purely 
analytical  formulation.  Diffraction  from  a  strip  grating  with  no  errors  has  been 
investigated  by  many  authors.  Most  of  the  analytical  solutions  have  involved  use  of  the 
Wiener-Hopf  method.  As  outlined  in  Appendix  F,  several  problems  exist  which  make  the 
use  of  the  Wiener-Hopf  method  inapplicable  to  the  problem  of  a  grating  with  errors  in  slit 
and  strip  widths.  Nevertheless,  a  concentrated  effort  to  develop  an  analytical  formulation, 
whether  based  upon  the  Wiener-Hopf  method  or  some  other  method  would  be 
worthwhile. 
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Appendix  A.  Variance  of  Linear  Array  of  Point  Sources 


This  appendix  derives  an  analytical  expression  for  the  variance  of  the  average 
power  pattern  of  a  hnear  array  of  point  sources.  To  keep  the  number  of  terms  at  a 
manageable  level  and  simplify  the  expressions,  the  variance  is  calculated  for  an  array  with 
an  infinite  number  of  point  sources. 

As  given  in  Equation  (2-2),  the  array  factor  for  the  array  of  point  sources  is  given 
by: 


N 

f(u)  =  ^ej"‘^" 

n=0 


(A-l) 


The  variance  of  the  average  power  pattern  is  given  by: 


G 
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The  second  term  is  the  square  of  the  average  power  pattern,  which  for  an  infinite  number 
of  slits/  strips,  is  given  in  Equation  (2-24): 
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where  O  is  the  characteristic  function  and 
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The  first  term  in  Equation  (A-2)  is  given  by: 
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The  quadruple  summation  can  be  expressed  as: 
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where 


Si  =  N  +  l  [for  indices  p  =  1  =  m  =  n]  (A-7) 

S2=N(N  +  1)  [for  indices  (p  “  1)  (m  =  n)]  (A-8) 

The  notation  (p  =  1)  (m  =  n)  refers  to  indices  such  that  p  =  1  and  m=  n  but  p  ^  mor  n 
and  1 mor  n.  Continuing  with  the  summations: 
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The  notation  (p  =  1)  m  n  refers  to  indices  such  that  p  =  1 ,  m  n,  and  p  mor  n. 
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Examining  the  above  equations,  the  following  are  true: 
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The  double  summations  can  be  expressed  as: 
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For  the  case  of  N  00,  the  following  is  true: 
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and  thus  there  are  no  contributions  from  S12,  S13,  S14,  and  S15.  A  similar  argument 
shows  that: 
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Also  for  the  infinite  case: 
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Now  consider  S4.  It  can  be  broken  into  24  cases  as  follows: 

Case  l:p>l>m>n  Case  13: 1  > p > n > m 
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(A-28) 

(A-29) 

(A-30) 

(A-31) 


Case  2:  p  >  1  >  n  >  m 
Case  3:  p  >  m  >  1  >  n 
Case  4:  p  >  m  >  n  >  1 
Case  5:p>n>l>m 
Case  6:  p  >  n  >  m  >  1 
Case  7:  m  >  1  >  p  >  n 
Case  8:  m  >  1  >  n  >  p 
Case  9:  m  >  p  >  1  >  n 
Case  10:  m  >  p  >  n  >  1 
Case  ll:m>n>l>p 
Case  12:  m  >  n  >  p  >  1 


Case  14: 1  >  p  >  m  >  n 
Case  15: 1  >  n  >  p  >  m 
Case  16: 1  >  n  >  m  >  p 
Case  17: 1  >  m  >  p  >  n 
Case  18:  l>m>n>p 
Case  19:  n>p>l>m 
Case  20:  n  >  p  >  m  >  1 
Case  21:  n  >  1  >  p  >  m 
Case  22:  n  >  1  >  m  >  p 
Case  23:  n  >  m  >  p  >  1 
Case  24:  n  >  m  >  1  >  p 


First,  consider  case  1,  p  >  1  >  m  >  n.  The  following  relations  are  true: 
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The  summation  becomes: 
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The  notation  S4  j  signifies  case  1  of  summation  4.  Next  consider  case  2,  p  >  1  >  n  >  m. 
For  this  case,  the  following  are  true: 


X„  =  Xo  +  L,+L2+-+L„  (A-38) 

X„  =Xo+L,+L2+-+L„+L„,i+-L„  (A-39) 

X,  =  Xq  +  Lj  +  Lj  +•  •  •  +L„  +  L„+i  +•  •  •  Ln  +  Ln+,  +•  •  •  h-L,  (A-40) 
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and 
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Now  consider  case  3,  p  >  m  >  1  >  n.  For  this  case,  the  following  are  true; 


-  Xq  +  Lj  +L2+---+Ln 

Xi  =Xo  +  Li  +L2+"-+Ln  +Ln+iH---Lj 

Xm  =^0  +L,  +L2+-"+Ln  +Ln+i+"-Li  +Li+i+---+L^ 

Xp  =  Xq  +  Lj  +  L2  +•  •  •  +L11  +  +•  •  ■  L,  +  L,^.i  +•  •  •  +L^  +  L^+i  +•  •  •  Lp 


and 
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The  summation  becomes: 
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For  case  4,  p  >  m  >  n  >  1,  the  following  are  true: 

Xj  =  Xq  +  Lj  +  Li2  +■  ■  •  "I'Ll 

Xn  =  Xq  +  Lj  +L2H  t-Lj  +Li^]H---Ln 
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Xm  -  Xq  +  Li  +  Lj  +•  •  •  +Li  +  L,+1  +•  •  •  Li  +  L,+1  +•  •  •  +L^ 


Xp  -  Xq  +  Li  +  L2  +•  •  •  +Li  +  Li+i  +•  •  •  L„  +  L„^.i  +•  •  •  +L„,  +  +•  •  •  Lp  (A-53) 
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which  is  the  same  as  case  3,  i.e. 

84,4=84.3  (A-56) 

For  case  5,  p  >  n  >  1  >  m,  the  following  relations  are  true: 

X„  =  Xo  +  L,+L2+--+L„  (A-57) 

Xi  =  Xq  +  Lj  +  Lj  +•  •  •  +L„  +  +•  •  •  L,  (A-58) 
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which  is  the  same  as  case  2,  i.e. 


^4,5  -  ^4,2 


For  case  6,  p  >  n  >  m  >  1.  The  following  relations  are  true: 

Xj  =  Xq  +  Lj  + 1^2  ‘  ‘ 

Xm  ^Xq  +  Li+Lj+'-'+L,  +Li+i+-”L„ 

=  ^0  +  Ll  +  ^2  +■  ■  ■  +Ll  +  Ll+1  +■  ■  ■  Lm  +  Lm+1  +' '  *  +^0 

Xp  =  Xq  +Li  +L2+-"+Li  +  Li+]+-"La  +  L^+i+"-+L„  +L„+j+***L 


and 


X„-X,+X„-X  =L, 


n+r 


•+Lp+L,+1+- 


•+U 


p-n  terms  tt"!  terms 


(A-62) 

(A-63) 

(A-64) 

(A-65) 

(A-66) 
p  (A-67) 

(A-68) 


A-10 


The  summation  becomes: 


S 


4,6 


N  p-1  n-i  m-1 


(p-n)+(m-l) 


p=3  n-2  m-1  1=0 


N  p-1  1-1  m-1  {  A  (  \ 

p=3  1-2  m-1  n=0 


(A-69) 


which  is  the  same  as  case  1,  i.e. 

84.6=84.1  (A-70) 

By  symmetry,  the  cases  when  "m"  is  the  highest  index,  i.e.  cases  7  -  12,  give  identical 
results.  Similarly,  the  cases  when  and  "n"  are  the  highest  indices,  i.e.  cases  13  -  24, 
result  in  the  complex  conjugates  of  the  above  cases.  The  final  result  is; 


84=  8  Re-1 


p=3  1-2  m-1  n=0 

(A-71) 


Let 


where 


and 


S4  —  8  Re{S43  +  S415  +84^} 

(A-72) 

N  p-1  1-1  m-1 

(A-73) 

p=3  1-2  m-1  n=0 

N  p-1  1-1  m-1  .  s 

=11X1®*“)'"  [«>■(“>] 

(A-74) 

p=3  1-2  m-1  n=0 


A-11 


(A-75) 


N  P  -1  1-1  m-1 


p=3  1-2  m-1  n=0 


The  following  relations  will  be  utilized  in  examining  these  summations: 


Let 


m=2 


n-1  . 

-V  ^ 

1- 

1/ 

/t" 

1 

~  1- 

1/ 

'“T 

-t^ 

m=0 

/t 

N 

N-l 

t(i- 

M 

3 

II 

M 

3 

II 

1- 

-t 

m=l 

m=0 

.  <(1- 

t^ 

(i- 

m _  ^ 

t= 

1- 

-t 

i-t 

0(u)  =  A 
a>(2u)  =  B 


Returning  to  S4,  843  can  be  written  as: 


where 


N  p~l  1-1  m-1  N  p~l  1-1 

S4.  =Ia'£a-'5;a"'£a-"  =Xa'Ia-‘Xa"si 

p=3  1=2  m=l  n=0  p=3  1=2  m=l 

N  p-1  N 

=  ^  A“'  s2  =  ^  A^  s3 

p=3  1=2  p=3 


m-1 


n=0 


(A-76) 


(A-77) 


(A-78) 


(A-79) 


(A-80) 


(A-81) 


A-12 


(A-82) 


1-1  1-1 


A  A' 

- rr +  T - 


1-A  (l-Af  (l-A)' 


A^p  riY.  Ap  2A'  ^  fM"'" 

(I-A)HaJ  (1-A)^  (1-A)’  (l-A)HAj 


(A-83) 


2  N 


2  N 


^43  “ 


(1-A)^ 


2  N 


(l-A)^ 


N 

^p-l-2  + 

P=1 

.2  f  N 


(1-Af^  (1-Ar^^ 


£pAi’-A-: 


p  J  Y  A' - 1  -  aH - 
(l-A)’  S 


(l-Af 

.  .,]  2A^(N-3) 

I  0-A)’ 


(A-84) 


A^  [KN  +  I)  1  A  rA  +  NA^^^^^-(N  +  l)A^^^^^  ^  ^^2 
^^~{l-An  2  yH-AY  {l-Af 


2A^  A-A^^"^^  ,  , 

+ - 5- - 1-A 

(1-A)^  1-A 


(i-ArL  y 

,  ,al  2A^(N-^ 


(A-85) 


(1-A)^ 


For  the  infinite  case. 


A-13 


(A-86) 


r  S4,  _ 

n!^(N+i)2  2(1 -A)^ 


Next,  S4U  can  be  expressed  as: 


N  p-1  1-1  m-l 


N  p-1  1-1 


p=3  1=2  m=l  n=0  p=3  1=2  m=l 

N  p-1  N 

=  ^AP^A-'s2  =  5^AP  s3 

p=3  1=2  p=3 


(A-87) 


where 


(A-88) 


1-1  *  ♦  y  V  A  *  r  1-1  1-1 

m=l  -1  ^  i  ^  Lni=0  m=0 


A*1  A*  A*^  ^ 
l-A’  (l-A-f'll-A-y 


(A-89) 


*  p-1 


A  *  p-1  A  *  P”^  A  * 

A-  l1  A-  ^  l],A-(p-2) 

1-A’ilr’  a|  (i-A-ns  aj  (i-A‘f 


AA*p  r  1  Y  .  A*p 


ll-ApUJ  |l-Ap  ll-Af(l-A)  ll-Af 


AA*(2  -  A  -  A*)  n  Y  A*A  ^(a*  -  a) 


|l-Ar 


(A-90) 


A-14 


S4b  “ 


A  A  *  N  A  *  N 

AA*(2-A-A*)  N 


2AA*  ^  A  A  ^(a  a) 
ll-Ap(l-A)^  |l-Ar  ^ 


(A91) 


which  becomes: 


S4b  = 


11 -Af 


Ip-1-2 


|l-Ap 


^pAP-A-2A^ 

P=1 


2AA*  A*A  ^(a*— a)  . 

+  - ^ - 1  JVap-i-a^ 

ll-Ap(l-A)  |1-A|"  1^ 


AA*(2-A-A*)(N-3) 

TaT 

(A-92) 


S4b  — 


AA* 

■N(N  +  1)  / 

A* 

|1-A 

2 

2 

1-A 

2 

L  2  J  |1-A|"[  (1-A)^ 

2AA*  A*A-'(A*-A)lrA-A('^"« 


ll-Ap(l-A)  |1-A|' 


-1-A"  - 


-A-2A"J 

AA*(2-A-A*)(N-3) 


ii-Ar 


(A-93) 


For  the  infinite  case, 


*4b  _ 

.  i\2 


N^-(N  +  ir  2|1-A|" 


(A-94) 


A-15 


Referring  to  Equations  (A-84)  and  (A-91),  the  only  contribution  as  N  — » oo  is  from  terms 
involving  ^p.  Next,  can  be  expressed  as: 

p 


N 

s2  =  ^APs3 


p=3 


(A-95) 


where 

m-1  A 

n=0 


(A-96) 


s2  = 


A 


1-A 


1-1 

IB- 


,m=0 


A 


1-A 


B 


1-B 


(b-'-i)- 


B/ 


1-B/ 


(A-97) 


which  becomes: 


A-16 


(A-lOO) 


s3  =  . 


A 

'(%) 

2  1 

[%) 

(1-A) 

■%: 

1  [ 

B/ 


B 


AB 


(1-A)(l-B)(l-V 


l-%  1-B 


1/V-  I 


B(p-2) 


A'  r(i-A)(i-9 


and 


s4  =  - 


A 

®/a 

B 

^  N  N  ' 

fB/fYAP-YB4 

(l-A)(l-%) 

1-B  ' 

v/A/  2^^ 

P=3  P=3  J 

_ AB 


N 


N 


/ArSA-ii 

P=3  P=3  J 


B 


N 


XpA'’-2£a'' 

P=3 


(A-101) 


While  these  summations  can  be  calculated,  the  above  expression  does  not  contain  a  term 
of  the  form  ^p  and: 

p 


Thus, 


lim-^4c.^  =  0  (A-102) 

(N  +  1) 


lim  lim  = 

n^-(N  +  1)  (N  +  1) 


4A^ 

(l-Af 


4AA* 


(A-103) 


Now  consider  S^).  It  can  be  broken  into  six  cases  as  follows: 

Case  4:  p  >  m  >1 
Case  5:  m  >  p  >  1 
Case  6:  m  >  1  >  p 


Case  l:p>l>m 
Case  2:  l>p>m 
Case  3: 1  >  m  >  p 


A-17 


First,  consider  case  1,  p  >  1  >  m.  The  following  relations  are  true: 


2Xp=2(Xo+Li+-+L„)  +  2(L^^,+-+L,)  +  2(L,^i+-+Lp) 


Xi  —  Xq  +  Lj  h  I-L^  +  H  Lj 


=  ^0  +  Lj  +  L2+-  •  •  +L^ 


and 


2Xp  -  X,  -  X„  =  2(l,^i  +•  •  •  +Lp) +(L„^i  +•  •  •  +L, ) 


p-1  terms 


1-m  terms 


The  summation  becomes: 


N  p-1  1-1 
p=2  1=1  m=0 


N 

sl  =  ^BPs2 
P=2 


where 


s1  =  £a-"=-^(a-‘-i) 

m=0 


l-A 


1=1 


and 


(A- 104) 

(A- 105) 

(A-106) 

(A-107) 

(A-108) 

(A-109) 

(A-110) 


A-18 


(A-111) 


or 


^10,1  “ 


1-A 


Sbp 


p=2 


B 


1-B 


(b-P-i) 


1 


JL 


which  becomes 


AB(N-I)  AB^(i-B^  Of  B_  A  ^ 
^^°’^“(1-A)(1-B)’*’  (1-A)(1-B)  [a-B  1-Ay 

A^b(i-A^“0 

"(i-a)"(a-b) 


and 


lim 

N-><» 


(N  +  lf 


=  0 


Now,  consider  case  2, 1  >  p  >  m.  The  following  relations  are  true: 

2Xp  =  2(Xo  +  L,  +• .  •  +L„)  +  2(l^,i  +•  •  •  +Lp) 

Xi  =  Xo  +  L,+---+L^  +  L^+i+-'-Lp  +Lp+,+"-L, 

Xm  =  ^0  +  L,  +  ^2+- "+Lni 


(A-112) 

(A-113) 

(A-114) 

(A-115) 

(A-116) 

(A-117) 


A-19 


and 


2Xp  -  X,  -  X„  =  L„,  +•  •  •  Lp  -  (Lp„  +.  •  ■  +L, ) 

> - ^  V - ^ - / 

p-m  terms  1-p  terms 


The  summation  becomes: 


N  1-1  P-1  ,  s  N  A  -\PPd 

*(l-p)  A  ' 


l=2p=lm=0  1=2  p=l^^ 


m=0 


1-1 


=  lATl^)  ^>  =  1^*  ^ 

1=2  p=l 


1=2 


where 


m=0 


(A- 11 8) 


(A-119) 


(A-120) 


(A-121) 


(A-122) 


A-20 


which  becomes 


aa*(n-i) 

AA’-^il 

1 

-A*‘  j 

/  1  _  1  1 

(1-aXi-a-) 

1  (1-a: 

U-A*  1-AV 

a’a-(i-a''-') 

(l-Af(A-A‘) 


and 


lim 

N  “—>00 


^10,2 

(N  +  ir 


=  0 


Next  consider  case  3, 1  >  m  >  k.  The  following  relations  are  true: 

2Xp=2(Xo  +  L,+-+Lp) 

Xj  =Xo  +  Li+-”+Lp  +Lp^.,+---Lm  +  L^+i+-"Li 
Xm  =  Xo  +  Li+L2+-"+Lp  +  Lp+i+---+L„ 

and 

2Xp-X,  -X„  =-  2(Lp,,+-L„)-(L„„+  -+L,) 

\ _ ^ — V - ^ 

m-p  terms  terms 


(A-123) 

(A-124) 

(A-125) 

(A-126) 

(A-127) 

(A- 128) 

(A-129) 


A-21 


The  summation  becomes: 


N  1-1  ra-1  \  .  N 


1=2  m=l  p=0 


,  1-1  /  o» 

lA'  S  It  Ib 

1=2  m=lV^  J  p=0 


1=2  m=lV^  / 


=  s1  =  £a*'s2 

^  1=2 


where 


si=yB--'’=^(B*-"-i) 
^  l-B 


1-1  fj3*Y'  n*  /  \  R*  1  £4 

^2  =  1  |t  S 

A  J  l-a  1  O  ni=A^  ''  1=1 


vin  m-1/  g» 


1=1  J 


B* 

l-B* 


A* 

1-A* 


(a 


and 


S  — ^ 


B 

"l-B* 


-1 


or 


SlA'J  —  • 


B*  A* 


l-B*  1-A* 


(N-.)-^(l-A-) 


B*  /B* 


l-B*  i_A*/ 

^  /b* 


B* 


l-B' 


1-A* 


(A-130) 


(A-131) 


(A-132) 


(A-133) 


(A- 134) 


A-22 


which  becomes 


A*B*(N-1)  B*A*^(i-A* _ )f  A*  A*  ^ 

Sio.3-^j_^*j^l_B*)+  (i_a*)(i-B*)  1,B*-A*  1-A*; 

b*V(i-an-0 

(i-b*)^(b*-a*) 


and 


Thus 


lim  =  0  (A-136) 

N^oo  (N  +  lf 


lim  T^^  =  0  (A-137) 

n^-(N  +  1) 


and  there  are  no  contributions  from  Sjg.  Since  Sj  j  is  the  complex  conjugate  of  S jg? 


lim  7-^T  =  0  (^-138) 

N-.00  (N  +  lf 


Collecting  the  non-zero  terms, 


(N  +  l)= 


!{if(»)r}= 


2  +  4Re- 


$2 


- --I- 4  Re - TT - r-rOI\.C - 

(l-<Df  (1-<D)(1-0*)  1-^> 


+  SRe- 


O 


(A-139) 


The  variance  is  given  by: 


A-23 


=  2  +  4Re 


+  4Re 


OO* 


+  8  Re - 

1-0) 


l  +  4Re 


I— I 

li-oj 


+  4 


Re^ 


O 


-]2 


1-0). 


or 


~  0)^  0)0)*  0) 

=l  +  4Re- - ^  +  4ReT - ^ - ;;y+4Re— — 

(l-O))'  (l-0))(l-0)*)  1-^ 


-4 


Re^ 


0) 


1-0). 


Simplifying, 


=l  +  4Re 


-0)^0)*  +0)0)*  -0)^0)*  +0)-0)'‘  -0)0)*  +0)^0) 


(l-4.)"(l-4>‘) 


-4^ 


Re< 


0) 


=l+4Re 


1-0. 
^(l-OO)*) 


-4 


MM 


Now, 


4Re 


0)(l-0)0)*)  0)(l-0)0)*)(l-0)*)  0)(l-0)0)*)-0)0)*+0)0) 


_ _  4  Rg _ 

(l-O)f(l-O)*)  (l-0))^(l-0)*)' 


=  4Re- 


(l-O)f(l-O)*) 


=  4Re 


^(O) +0)*)(l  -  0)0)*)  -  0)0*  +  0)0)* 


(l_0)f(l_0)*) 

20)  +  20)*  -  20)^0)*  -  20)0)*^  -  40)0)*  +  40^0)* 
(l-0)f(l-O*)' 


(A- 140) 


(A- 141) 


(A-142) 


(A-143) 


A-24 


and 


„  <D  „  a>(i-o*)  i(o+a)*)-<DO* 

" ^^(i-<i>)(i-o*)  ~  (i-a>)(i-a)*) 


which  leads  to: 


Re 


O 


1-<D 


=  4 


i(3)  +  <E)*)-0<D’ 
(l-4>)(l-4«-) 


^2  +^*^20$*  -40^0*  -400*^  +  4<[>^0*^ 


Substituting  Equations  (A- 143)  and  (A- 145)  into  Equation  (A- 142)  yields: 


a"=l  + 


2<I)  +  2<1)*  +20^0*  +20a>*^  -bOO* 


(l-4>)"(l-<I>-) 

(1_0)^(i-4)")^  +  2i1i+24>‘+24>’®‘  +  24>®"^-6**"-4>^ 


or 


a 


2 


l-20C&*+0^0*^  (l-0$  ) 

(l  -  o)^  (l  -  <D*f  (l  -  0)(l  -  O*)  ^ 


(A- 144) 


(A-145) 


2  (A-146) 


(A-147) 


A-25 


Thus  the  variance  is  equal  to  the  square  of  the  average  power  pattern  and  the  standard 
deviation  is  equal  to  the  average  power  pattern,  i.e. 


Equations  (A- 147)  and  (A- 148)  are  valid  for  u  0.  For  the  case  of  u  =  0, 


N  N  N  N 
p=0  1=0  m=0n=0  J 


=--^e{(n+i)‘}=(n+i)’ 

(n+i) 


and 


=  (n+i)^-(n+i)^ 

=  0  ,  for  u  =  0 


(A- 148) 


(A-149) 


(A-150) 


(A-151) 


A-26 


Appendix  B.  Modal  Formulation  for  Infinite  Grating  with  No 

Errors 


Conversion  of  Integral  Equation  to  Modal  Equation 

Setting  the  origin  at  the  mid  point  of  the  first  slit  or  strip,  i.e.  1/2  (Xq  +  Tq)  =  0, 

the  scattered  fields  for  a  grating  with  no  errors  are  given  by  (Equation  (3-51)): 


S  =  -X j".  Is (X")  (k^(x-x"-nT)^-^y^)dx" 

^  n=0  2 

=  “Sjtls(x")  —  Hf\k^(x-x"-nT)'+y")| 


n=0  2 


d(odx" 

(B-1) 


where 


and 


(?'{f(x)}  =  Jf(x)e  ■'^’'dx 
f(x)  =  {f(x)}e^“'‘d(0 


(B-2) 

(B-3) 


are  the  Fourier  and  inverse  Fourier  transform  relationships.  Using  the  Fourier  Transform 
(2:58): 


(B-4) 


the  scattered  fields  become: 


B-1 


3  = 

^  n=0  2 


JCOX 


JnkTsin6i 


^/i?^ 


dcodx' 


(B-5) 


k  rf  _  -ejco(x-x")g-j|y|V?=;? 


2k 


Vk^  -0)^  „=o 


^ej"kTs‘n0ie-j“nT  dcodx' 


or 


S  =  — Is(x")  f 

2k  j 


Ja)(x-x")g-j|y|Vk^-w^  27t 

Vk^-w^  ~ 


^  5(o-(ksin6i  +-2^))do)dx' 


V  jco(x-x")  -j|y|Vk^-<o'  - 

=  J  ■^Xs(<->-PJdo.dx- 

A  •’"T  '  Vk  -0)  n=-oo 


(B-6) 


where 


o  ,  •  n  27m 

Pn=ksmei+— 


(B-7) 


Using  the  sifting  property  of  the  delta  function,  the  scatterred  fields  become: 


;=iV 


"  -W  jPn(x-X")  -jlylVk^-Pn^ 

Is(x") - .  - dx' 


T  ^ 


n=-oo 


Vk^ 


1,  ~  *w  ^jPn(x-x")p-jiy|pn 

=1 1  ri.(x")‘  p/  dx- 


(B-8) 


where 


Pn=>/k^-Pn 


(B-9) 


Equation  (B-8)  can  be  rewritten  as: 


B-2 


fory>0 


(B-10) 


where 


S  = 


E®"' 

n 

E®«' 


.XPnX-Pny) 


j(Pn*+Pny) 


L  n 


fory  <0 


(B-11) 


Equations  (B-10)  and  (B-11)  are  equivalent  to  modal  formulations  found  in  many 
references  ((11:189-194)  for  example).  For  the  TEz  mode,  coefficients  of  reflection  (R) 
and  transmission  (F)  are  given  by: 


R  =  l-B 
r  =  Bo 


0 


(B-12) 


and  for  the  TMz  mode: 


R  =  -Bo 
r=i-Bo 


(B-13) 


Born  Approximation 

Using  the  Bom  approximation  as  given  in  Equation  (4-9): 


I,(x")  =  cosei 


(B-14) 


then, 


B-3 


B 


n 


kcOsGj  1*2  gXksinei-P„)x" 

n  T  j-f 


P,T 

kcosSi 


PnT 


^  e‘'T"  dx" 

w 


wkcosGj  sin(-°^) 

■T’  RTiSL 

Kn^  T 


with 


Bo 


w 


T 


Edge  Condition  Approximation 

Using  the  edge  condition  approximation  as  given  in  Equation  (5-3): 


then 


(B-15) 


(B-16) 


(B-17) 


(B-18) 


B-4 


To  determine  the  value  for  the  constant  "C",  use  the  integral  equation  relating  the 
boundary  conditions: 


^  ^  jkxsinOi 

k 


(B-19) 


Multiply  both  sides  by  Ig*  and  integrate  over  -w/2  to  w/2  to  obtain: 


or 


dx 


Ck7iv'Jo(^)  f' 


f  (w^f-(x)^ 


T  „  Pn  J-f[(W/2)2_(x'f]^^ 


n  r-ji 


(B-20) 


C7l  = 


n  Hn 


C  = 


ktt 


1-1 


T  "p 

n  ^^n 


(B-21) 


B-5 


Appendix  C.  Simplification  of  Integral  Used  in  Edge  Condition 

Approximation 


This  appendix  simplifies  the  following  integral  used  in  the  edge  condition  approximation: 
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Appendix  D.  Variance  of  the  Average  Power  Pattern  of  a  Strip 
Grating  -  Born  Approximation 


This  appendix  derives  an  analytical  expression  for  the  variance  of  the  average 
power  pattern  of  a  strip  grating  computed  using  the  Bom  approximation.  As  in  the  case 
of  the  array  of  point  sources,  the  variance  is  calculated  for  a  grating  with  an  infinite 
number  of  slits/  strips  to  keep  the  number  of  terms  at  a  manageable  level  and  simplify  the 
expressions. 

I 

The  variance  of  the  average  power  pattern  is  given  by: 


=  (^;^^e{is(u)|"}-[|s(u)P 


(D-l) 


The  second  term  is  the  square  of  the  average  power  pattern,  which  for  an  infinite  number 
of  slits/  strips  and  using  the  Bom  approximation ,  is  given  in  Equation  (4-35): 
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The  first  term  in  Equation  (D-l)  is  given  by: 
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Following  the  same  procedure  used  in  Appendix  A  for  the  array  of  point  sources,  this 
becomes  for  the  case  of  N  ^  oo : 
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Substituting  Equations  (D-2)  and  (D-4)  into  Equation  (D-1),  the  variance  becomes: 
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To  simplify,  let: 


a  =  E 
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The  second  and  third  terms  of  Equation  (D-5)  become: 
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The  fifth  term  of  Equation  (D-5)  can  be  expressed  as: 
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Substituting  Equations  (D-7)  and  (D-8)  into  Equation  (D-5),  the  variance  becomes: 
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or  after  substituting  Equation  (D-8)  for  the  last  term: 
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and  the  standard  deviation  is  equal  to  the  mean,  i.e. 
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The  above  equations  are  valid  for  u  0.  Substituting  u  =  0  into  Equation  (D-3)  yields: 
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Sustitution  of  the  last  two  equations  into  Equation  (D-1)  yields: 
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which  goes  to  infinity  as  N  ->  <».  On  the  other  hand: 
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Appendix  E.  Variance  of  the  Average  Power  Pattern  of  a  Strip 
Grating  -  Edge  Condition  Approximation 


This  appendix  derives  an  anal5rtical  expression  for  the  variance  of  the  average 
power  pattern  of  a  strip  grating  computed  using  the  Edge  condition  approximation.  As  in 
the  cases  of  the  array  of  point  sources  and  the  grating  using  the  Bom  approximation,  the 
variance  is  calculated  for  a  grating  with  an  infinite  number  of  slits/  strips  to  keep  the 
number  of  terms  at  a  manageable  level  and  simplify  the  expressions. 

As  given  in  Equation  (D-l)The  variance  of  the  average  power  pattern  is  given  by; 
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The  second  term  is  the  square  of  the  average  power  pattern,  which  for  an  infinite  number 
of  slits/  strips  and  using  the  edge  condition  approximation ,  is  given  in  Equation  (5-32): 
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The  first  term  in  Equation  (E-1)  is  given  by: 
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Following  the  same  procedures  used  in  Appendix  A  for  the  array  of  point  sources  and 
Appendix  D  for  the  grating  using  the  Bom  approximation,  this  becomes  for  the  case  of 
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Substituting  Equations  (E-2)  and  (E-4)  into  Equation  (E-1),  the  variance  becomes: 
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To  simplify,  let: 
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The  second  and  third  terms  of  Equation  (E-5)  become,  using  Equation  (D-7): 
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The  fifth  term  of  Equation  (E-5)  can  be  expressed  as,  using  Equation  (D-8): 
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Substituting  Equations  (E-7)  and  (E-8)  into  Equation  (E-5),  the  variance  becomes: 
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or  after  substituting  Equation  (E-8)  for  the  last  term: 
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and  the  standard  deviation  is  equal  to  the  mean,  i.e. 

E-3 


The  above  equations  are  valid  for  u  0.  Substituting  u  =  0  into  Equation  (E-3)  yields; 
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and  from  Equation  (5-30): 
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Substituting  the  last  two  equations  into  Equation  (E-1)  yields: 


Appendix  F.  Problems  with  Finding  an  Analytical  Solution 


As  stated  in  Chapter  I,  the  problem  of  diffraction  by  a  strip  grating  has  been 
investigated  by  many  authors.  The  pioneering  works  in  solving  this  problem  analytically 
was  accomplished  by  Baldwin  and  Heins  (5)  and  Weinstein  (29)  using  the  Wiener-Hopf 
method.  Several  problems  exist  which  make  use  of  the  Wiener-Hopf  approach 
inapplicable  to  the  problem  of  a  strip  grating  with  random  errors  in  strip  width  and 
spacing.  First,  the  mathematical  formulation  of  the  problem  of  diffraction  by  a  grating 
containing  an  infinite  number  of  strips  and  slits  results  in  an  infinite  set  of  integral 
equations.  The  periodic  property  of  the  solution  allows  the  infinite  set  to  be  replaced  by  a 
single  integral  equation  valid  over  one  period.  Once  the  solution  over  one  period  is 
obtained,  the  solution  to  other  periods  is  obtained  through  phase  adjustment  (Floquet's 
theorem).  Second,  an  analytical  solution  using  the  Wiener-Hopf  method  is  possible  only 
when  the  width  of  the  strips  is  equal  to  the  width  of  the  apertures.  As  a  result,  use  of  the 
Wiener-Hopf  approach  is  inapplicable  to  the  problem  of  finding  an  analytical  solution  to 
the  problem  of  diffraction  from  a  grating  with  random  strip  and  slit  widths. 
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